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The  aim  of  this  project  was  to  obtain  a  more  complete  undCTStanding  of  the  physical  phenomena 
occuning  in  piezoelectric  transducers,  with  the  aim  of  obtaining  optimal  designs  and  improved 
manufacturability.  Of  key  importance  in  design  are  conversion  efficiency,  bandwidth,  and  control  of 
spurious  modes.  The  transducer  structures  considered  are  1:3, 0:3, 2:2,  and  1:3:1.  Contributions  to 
the  advancement  of  1 :3  modeling  and  fabrication  are:  (i)  finite  element  modeling  of  composite  lattices 
with  complex  unit  cells,  (ii)  transmission  line  and  reciproci^  modeling  of  distributed  jxmod 
composites,  (iii)  fabrication  of  test  transducers  for  verification  of  the  modeling  principles.  New 
naodeUng  theory  for  0:3  composites  aims  at  replacing  the  empirical  “cubes”  model  with  a  nnodel 
consisting  of  nearly  close-packed  spherical  particles.  This  approach  estimates  mrae  closely  the 
performance  of  0:3  composites  and  points  the  way  for  further  studies.  The  2:2  structures  coiisidered 
are,  in  fact,  multilayered  electroded  transducers  with  a  single  piezoelectric  material.  But  multiple 
piezoelectrics  and  dtemating  piezo  and  nonpiezo  layers  may  also  be  considered.  A  procedure  has  been 
(tevelofxxl,  and  confirmed  experimentally,  for  tailtning  the  transducer  time-domain  response  by 
af^lying  multiple  excitations.  Hnally,  1:3:1  transducers  with  a  periodic  additicm  of  empty  voids  have 
b^  developed.  This  version  produces  a  more  uniform  thickness  mode  displacement  and  improved 
coupling,  pr^cted  finom  finite  element  modeling  and  confirmed  experimentally. 
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1.  INTRODUCTION 

The  objective  of  this  project  was  to  obtain  a  more  complete  understanding  of  the 
physical  phenomena  occurring  in  piezoelectric  composite  transducers.  Techniques  for 
controlling  parasitic  resonance  mode  activity  in  1:3  composite  transducers  have  been 
developed.  The  mecheinisms  governing  the  behavior  of  0:3  transducers  have  been  considered 
in  detail.  The  concept  of  using  transducers  with  separately  excited,  but  mechanically 
coupled,  layers  has  been  proposed.  Finally,  an  enhanced  1:3  composite  transducer,  with  a 
compressible  third  phase,  has  been  developed. 

2.  1:3  COMPOSITE  TRANSDUCERS 

This  type  of  composite,  comprising  piezoceramic  pillars  embedded  in  a  polymer 
matrix,  has  applications  in  both  active  and  passive  modes  of  operation.  The  present 
work  has  concentrated  on  identifying  the  influence  of  pillar  geometry  and  configuration  on 
device  performance. 

2.1.  Floquet  Analysis 

In  this  approach  the  material  parameters  (density  and  stiffness)  are  written  as 
arbitrary  periodic  fimctions  of  the  spatial  coordinates.  In  a  1-3  composite  the  material 
properties  are  periodic  in  a  plane  (xy)  and  uniform  along  the  perpendicular  axis. 

The  advantage  of  this  approach  is  that  it  permits  a  full  treatment  of  an  infinite  1-3 
composite.  Simple  approximations  for  small  volume  fraction  (<  10%)  and  high  volume 
fraction  (>  90%)  of  ceramic  give  excellent  physical  insight  into  the  behavior  of  waves  and 
vibrations  in  these  materials,  and  facilitate  the  development  of  design  concepts  based  on 
physic2d  principles.  Within  the  above  limitations  this  approach  permits  a  full  development 
of  wave  dispersion  curves  and  latersJ  standing  wave  resonances  at  all  of  the  material 
stopbands.  The  approach  is  not  practical  for  intermediate  volume  fractions,  where  the  finite 
element  method  (2-dimensional  periodicity)  and  the  wave  transmission  matrix  method 
(1 -dimensional  periodicity)  must  be  used  to  flesh  out  the  details.  Floquet  theory  can, 
in  principle,  be  applied  to  the  analysis  of  a  composite  plate,  but  this  is  not  practical 
because  of  the  complexity  of  the  equations.  Finite  element  emalysis  must  be  used  to  make 
a  connection  between  Floquet-based  concepts  and  the  behavior  of  1-3  composite  plates, 
especially  those  having  intermediate  volume  fractions  of  ceramic  in  their  composition. 

(a)  Basic  Wave  Equation 

Previous  analysis  and  experiment  show  that  lateral  resonances  in  1-3  composites  are 
basically  shear  waves  polarized  parallel  to  the  length  direction  of  the  ceramic  elements. 
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A  simple  periodic  wave  equation  can  be  developed  for  model  calculations  based  on  this 
approximation. 

Define  a  coordinate  system  with  x  and  y  axes  in  the  basal  plane  and  the  z  axis  parallel 
to  the  ceramic  elements.  The  displacement  field  of  the  z-polarized  vibration  is 


tir(x,y) 


Since  the  z  displacement  is  a  function  of  both  x  and  y,  the  strain  componenets  are 


Szy 


_  1  _idu, 

2  2  dx 

_  1  1  au, 

2  2  dy 


If  the  two  composite  phases  are  isotropic  the  stresses  are 


Txx  =  =  C44 


dx 


Tiy  =76=  C44 


The  dynamic  equation 


dTij  2 

-a-  =  -p- 


then  yields 

or 

/ d  d  d  d\  2 

Consider  a  composite  with  different  periods  along  the  horizontal  (i)  and  vertical  (y)  axes, 
as  shown  in  part  (a)  of  Fig.  1,  from  the  proposal. 
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(b) 
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saw 
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Figure  1 

Alternative  1-3  composite  lattices 
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In  this  case  C44  and  p  are  doubly  periodic  functions  of  x  and  y,  with  period  dx  along 
X  and  d,  along  y.  Write 


C44(x,  y)  =  C44  +  6  C44(x,  y) 


p{x,y)  =  /)  +  ^p(x,y) 


(2) 


where  the  barred  quantities  are  the  (constant)  spatial  average  values  and  6c44(x,  y),  6p{x,  y) 
are  doubly  {>eriodic  functions,  as  noted  above. 

Substitution  of  Eq.  (2)  into  Eq.  (1)  gives 


{cn^ly  +  poj^)  ttj  =  -6W{x,y)u: 


(3a) 


where 


6W{x,y)  =  C44(a:,y)V^  +ic44(x,y)V 


d  d  d  d 


(b)  General  Floquet  Solution 

The  general  solution  to  Eqs.  (3a)  and  (3b)  has  the  form 


«z(x,y)  =  -  r) 

FFi.n 


where 


r  =  XX  +  yy 


L 


(4a) 


(46) 


with  m,n  0,  ±1,±2,  ±3,  etc.  and  0oz,0Qy  arbitrary.  Substitution  into  Eq.  (3a)  yields 


'^{pu'^  -C440l,n)'i:  =  -^Wix,y)Uz  (5) 

m,n 
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The  problem  is  to  calculate  the  unknown  space  harmonic  amplitudes  Omn  in  £q.  (4a). 
In  Eq.  (5)  u,  on  the  LHS  is  written  as  in  Eq.  (4a),  and  on  the  RHS  it  is  written  with  m,  n 
changed  to  p,  q.  Note  that 


r-  ay  (e‘t. A- 

Jo  Jo  \  /  \  /  y  dxdy  r  =  m,s  =  n 


because  of  the  definition  of  in  Eq.  (4b).  Multiplying  both  sides  of  Eq.  (5)  with 


and  integrating  over  a  imit  cell  in  Fig.  1(a)  yields  a  single  term  on  the  LHS  and  gives 


^44  0ra)  ^rs  —  ^  ^ra,pqO’pq  (7) 

P9 

with 

^ ~  lo  lo  *  {^ra  ~  ipg)  '  ^ 

where  rs  nms  over  all  positive  and  negative  integers,  including  zero.  Equation  (7)  is  an 
infinite  set  of  linear  equations  defining  the  space  harmonic  eimplitudes  in  Eq.  (4). 

(c)  Coupled  Mode  Approximation 

Prom  Eq.  (4b)  fi^a  LHS  of  Eq.  (7)  is 

ffh  =  (^0.  +  +  (^0,  +  (8) 

The  equation  set  (7)  has  an  infinite  set  of  LHS's,  each  with  different  indices.  In  E^q.  (4a) 
the  leading  term  has  m  =  n  =  0,  and  the  corresponding  LHS  in  Eq.  (7)  is 


(^*  -  C44$lo)  flOO,  with  000  =  {0lx  +  0Oy) 


8 


For  a  spatially  uniform  medium  the  RHS’s  of  Eqs.  (5)  and  (7)  are  zero,  so  that  aoo 
is  the  only  nonzero  amplitude.  For  a  periodic  medium,  all  space  harmonics  exist,  but 
some  are  larger  than  others.  It  can  be  shown  that  when  two  of  the  LHS’s  of  Ekj.  (7)  axe 
simultaneously  close  to  zero,  the  corresponding  space  harmonic  amplitudes  are  large.  This 
corresponds  to  the  center  of  a  stopband.  If  the  RHS’s  of  Eq.  (7)  are  small,  all  of  the  other 
space  harmonic  amplitudes  are  small  and  can  be  neglected. 

In  the  above  situation  Eq.  (4a)  is  approximated  by  trvmcating  to  the  two  “near 
resonance”  space  harmonics — i.e.,  harmonics  for  which  the  LHS’s  axe  close  to  zero.  This  is 
called  the  coupled  wave  approximation.  In  this  case  Eq.  (7)  reduces  to  two  or  more  linear 
equations  in  the  same  number  of  unknowns.  For  example, 

(9) 

^440^k)  ~  ^jkfhi^hi 

The  characteristic  determinant  of  this  set  gives  the  dispersion  relation  near  the  hi,jk 
stopband.  That  is 

(^2  -  C44^L)  -  C440^ik)  =  KurJkKjkM  (10) 

Solution  of  Eq.  (10)  for  u;  versus  0o  in  Eq.  (4b)  defines  the  dispersion  cinrve  of  a 
stopband  centered  on  the  frequency  where  the  two  factors  of  the  LHS  of  Eq.  (10)  are 
simultaneously  zero.  This  is  a  characteristic  feature  of  Floquet  theory.  Stopbands  occm 
at  frequencies  when  two  or  more  of  the  f2w:tors  on  the  LHS  of  Eq.  (7)  axe  simultaneously 
zero. 

In  composite  theory  the  most  important  stopbands  relate  to  space  harmonics  that 
propagate  along  either  the  x-axis  or  the  y-axis — that  is,  with  either  m  or  n  equal  to  zero  in 
Eq.  (4b).  For  small  volume  fraction  composites,  where  the  coupled  mode  approximation  is 
valid,  it  can  be  seen  that  when  two  of  the  factors  on  the  LHS  of  Eq.  (7)  are  simultaneously 
zero  there  occurs  a  Bragg  scattering  (or  a  stemding  wave  resonance)  along  either  the  x 
or  the  y  axis.  The  resonances  occur  between  adjacent  rows  of  ceramic  elements  aligned 
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perpendicular  to  the  propagation  direction.  In  Fig.  1(a)  it  is  clear  from  the  dimensions  that 
x-directed  and  y-directed  resonances  do  not  occur  at  the  same  frequency.  However,  it  can 
be  shown  by  evaluating  the  coupling  constants  in  the  RHS  of  Eq.  (7)  that  these  standing 
waves  are  not  coupled.  The  resonances  are  therefore  uncoupled  one-dimensional  standing 
waves.  This  result  has  been  found  to  be  consistent  with  a  finite  element  analysis.  In  Fig. 
1(c)  X-  and  y-directed  standing  wave  resonances  now  exist  at  a  single  frequency.  But  now, 
by  contrast,  they  are  coupled,  so  that  the  resonances  axe  two-dimensional  standing  waves. 

2.2.  Finite  Element  Analysis 

It  has  proven  possible  to  obtain  a  match  between  analytic  (Floquet  etc.)  and  finite 
element  results.  A  10%  composite  was  analyzed  using  finite  element  analysis  and  the 
frequencies  of  the  thickness,  and  first  two  lateral  modes  were  determined  to  be  605  kHz, 
793  kHz,  and  1109  kHz,  respectively.  The  measured  frequencies  of  these  modes  were  621 
kHz,  802  kHz  and  1033  kHz.  Hence,  good  agreement  is  evident.  Furthermore,  the  finite 
element  displacement  plots,  shown  in  Fig.  2,  indicate  that  the  first  and  second  lateral 
modes  correspond  to  resonant  activity  principally  parallel  and  diagonal  to  the  planes  of 
pillars.  These  modes  are  better  illustrated  in  Fig.  3  which  presents  cross-sections  through 
the  thickness  direction  of  a  2%  composite.  The  first  mode  is  shown  along  a  section 
parallel  to  a  pl2ine  of  pillars.  The  second  mode  is  shown  for  a  diagonal  cross-section. 
In  accordance  with  the  modes  predicted  using  Floquet  theory,  the  first  mode  consists 
of  one  parallel  propagating  wavelength  and  the  second  mode  consists  of  two  diagonally 
propagating  wavelengths.  However,  when  larger  volume  fractions  or  wider  pillar  aspect 
ratios  are  considered,  the  complexity  of  the  resonant  behavior  necessitates  the  exclusive 
use  of  numerical  techniques. 

Some  novel  pillar  spacing  geometries  have  been  considered.  Figure  1(a)  illustrates  the 
lowest  frequency  bandedge  resonance  for  one  particular  configuration.  It  possesses  half 
wavelength  spacing  in  the  horizontal  direction  zmd  quarter  wavelength  spacing  adong  the 
vertical  direction.  Therefore,  alternate  vertical  columns  vibrate  zmtiphase  and  hence  the 
net  piezoelectric  coupling  is  zero.  A  similar  effect  occurs  at  the  next  higher  resonance, 
with  X  spacing  along  the  horizontal  axis  and  A/2  spacing  along  the  vertical  axis.  This 
structiure  was  analyzed  using  finite  element  analysis  and  results  obtmned  from  a  similar 
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Figure  2 

Modal  displacement  plots  of  a  10%  composite 


First  Mode — Parallel,  Wavelength  Resonance 


Second  Mode — Diagonal,  Double  Wavelength  Resonance 

Figure  3 

Cross-section  of  a  2%  Composite 
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regular  square  section  geometry  were  compared. 

Figxue  4  illustrates  the  fimdamental  thickness  mode,  an  extremely  weakly-coupled 
antisymmetric  mode  and  the  first  and  second  lateral  modes  corresponding  to  the  regular 
square  section  composite.  The  electromechanical  coupling  coefficients  associated  with 
each  of  these  modes  are  0.178,  0.0,  0.300  and  0.212,  respectively.  However,  when  the 
rectangular  configuration  was  analyzed,  the  “weakly”  coupled  mode  was  no  longer  weak 
and  had  split  into  two  separate  resonances  corresponding  to  length  and  width  dimensions, 
respectively.  The  higher  frequency  mode  remained.  All  these  mode  are  illustrated  in  Fig. 
5.  The  coupling  coefficients  associated  with  these  modes  are  0.121,  0.287,  0.216  and  0.210, 
respectively.  Consequently,  the  rectangular  geometry  has  produced  additional,  and  initially 
unexpected,  modes.  These  modes  are  sufficiently  strongly  coupled  electromechanically  that 
they  rule  against  the  use  of  this  design  in  practical  applications. 


2.3.  Wave  Transmission  Matrix  Technique  for  Distributed  Dimension  Composites 

Although  finite  element  analysis  provides  the  most  versatile  means  of  analysis  of 
composite  structures,  it  is  very  computationally  intensive.  Furthermore,  it  is  unsuited 
to  the  analysis  of  more  than  a  few  ceramic/polymer  cells.  If  transverse  wave  propagation 
is  of  primary  concern  and  the  approximation  made  in  using  a  one-dimensional  stack  of 
ceramic/polymer  plates  is  acceptable,  very  much  more  efficient  analysis  techniques  are 
available.  The  method  employed  here  was  developed  and  described  by  Wang  [1].  Consider 
the  layered  structure  illustrated  in  Fig.  6.  Using  a  knowledge  of  layer  dimensions  and 
zicoustic  characteristics,  it  is  possible  to  derive  a  matrix  describing  a„,  b„  in  terms  of  a„, 
bn-\ .  This  resvdts  in  the  following  expression  after  considerable  algebraic  manipulation, 
which  is  detailed  elsewhere  [2]. 


fln-i 

bn-l 


(A  B 
VC  D 


(11) 
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Thickness 


Weak  Mode 


First  Lateral  Second  Lateral 


Figure  4 

Modal  Displacement  Plots  of  a  Square  Section  0.5%  Composite 


Thickness 


Length  direction  lateral  mode 


Width  direction  lateral  mode 


Secondary  lateral  mode 


Figure  5 

Modal  Displacement  Plots  of  a  Rectangular  Section  0.5%  Composite 
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Figure  6 

One-dinMnsional  Layered  Medium 
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where 


_  Aj/j) 

titi 


A  =  - - C 

hh 


B  =  —  g«(*i^t+Aa<2)  ^  21  (f  2L22.^  g»(*i<i-A2<2) 

L<1<2  <1  V  <2  7 

2l  (t^  —  22-H^  «(Airi+ij/j) 

<2  V  <1 


C  =  g— •(Ai<i~Aj/j) 


£>  =  - 


J1  g-.(t.<i-Aa/a)  -  f<2  _  2121^  _  2122^  g-.(A,r, 

/l<2  V  /  V  ^2  / 


+*2^2) 


(12) 


Using  this  information,  it  is  possible  to  derive  the  Bloch  wave  dispersion  characteristics 
using  the  following  expression. 


K 


\^A  +  n) 


(13) 


Additionally,  the  transmission  coefficient  of  a  wave  transmitted  through  a  stack  containing 
numerous  layers  may  be  derived.  Assuming  that  there  is  no  reflection  from  the  last  layer 
(propagation  medium),  6„  =  0.  Hence  the  following  expression  is  obtained. 


Therefore 


Oo 


1 

A' 


(14) 


where  A'  corresponds  to  the  element  in  the  A  position  in  the  n***  power  matrix, 
evaluated  from  Eq.  (12)  for  the  specified  Bloch  wave  solution.  This  technique  permits 
the  lateral  transverse  wave  propagation  characteristic  to  be  evaluated  for  a  diverse  range 
of  ceramic/polymer  configurations. 
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The  most  promising  method  identified  for  the  suppression  of  undesirable  lateral 
periodic  rescmances  involves  the  use  of  diverse  dimensions  (distributed  periods)  between 
ceranoic  pillars.  In  diis  repent  recent  theoretical  and  experimental  work  is  reviewed.  Senne 
of  diis  wenk  has  recently  been  published  [3]. 

(a)  Results  Obtained  Using  Modeling  Techniques 

It  has  been  established  that  the  degree  of  attenuation  experienced  by  laterally 
propagating  waves  in  a  1:3  composite  structure  is  a  function  of  the  variance  in  the 
distribution  used  to  determine  pillar  spacings  [1].  The  attenuation  is  also  a  function  of 
volume  fraction  and  it  should  be  noted  that  variance  may  be  applied  to  a)  the  pillar  cross- 
sectional  dimensions,  b)  the  inter-pillar  spacing,  or  c)  both  a)  and  b)  simultaneously. 
Alternatively,  diversity  of  material  characteristics  may  be  introduced — e.g.  as  alternating 
slices  filled  with  hard  and  soft  polymers. 

Figure  7  illustrates  the  dependence  of  wave  attenuation  on  frequency  through  a  400 
layer  10%  one  dimensional  composite  compridng  layers  of  PZT-SH  and  HYSOL  RE  2039 
epoxy.  The  regular  device  consists  of  0.2  mm  ceramic  strips  with  1.8  mm  polymers  strips 
between  them.  It  is  immediately  sqiparent  that  increasing  variance  (applied  simultaneously 
to  bodi  ceramic  and  polymer  dimensions)  results  in  increasing  wave  attenuation.  Using  the 
same  basic  geometry,  a  comparison  is  made  between  the  effect  of  applying  a  25%  variance 
only  to  the  ceramic  and  then  only  to  the  polymer.  Notice  that,  if  the  composite  is 
manufactured  using  the  “slice  and  fill”  technique,  then  it  is  normally  simpler  to  apply 
variance  to  the  pillar  dimensions  since  the  saw  kerf  is  restricted  to  the  blade  width — 
assuming  single  saw  passes  with  one  bla^  size.  Figure  8  indicates  that  for  this  geometry  a 
far  greater  attenuation  effect  is  obtained  by  applying  variance  to  the  polymer.  Since  the 
pol3rmer  dimension  is  greater,  and  the  speed  of  sound  slower  in  the  polymer,  with  respect 
to  die  ceramic,  a  given  percentage  change  in  the  dimension  of  the  polymer  will  result  in  a 
greater  phase  shift. 

If  the  opposite  end  of  the  volume  fraction  range  is  now  considered,  different  results 
are  obtained.  Hgure  9  illustrates  the  relative  effects  of  25%  variances  aj^lied  separately  to 
the  ceramic  and  ptdyroer  phases.  In  this  case  applying  variance  to  the  ceramic  dimension 
has  a  greater  effect  since  the  wave  is  within  the  ceramic  phase  for  a  greater  propention  of 
the  time.  An  alternative  simple  procedure  for  disrupting  the  stopband/passband 
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Figure  7 

Lateral  tranamianon  through  a  400  layer  structure  of  alternating  PZT-5A  ceramic  and 
Hyaol  RE  2039  epoxy. 

Stngte  period  (0.2  mm  ceramic,  1.8  mm  po^mer) - 

5%  Variaaoe  in  dimenaiona  of  both  phasM  - 

25%  Variance  in  dimenaiona  of  both  phases  — . 

50%  Variance  in  dimenaions  of  both  phases - 
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Figure  8 

Lateral  transmisaon  through  a  400  layer,  10%,  structure.  Variance  applied  to  different 
phases 

25%  Variance  ai^lied  to  ceramic  dimensions  - 

25%  Variance  applied  to  p<^]nDer  dimensi<ms  . . 

25%  Variance  applied  to  dimensions  of  both - 

ceramic  and  polymer 
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Figure  9 

Lateral  transmission  through  a  400  layer,  90%,  structure.  Nominal  dimensions;  1.8  mm 
PZT-5A  ceramic  and  0.2  mm  Hysol  RE2039  epoxy. 

25%  Variance  applkd  to  ceramic  dimensicms  - 

25%  Variance  applied  to  polymer  dimensions  - - 

25%  Variance  applied  to  dimensions  ci  both - 

ceramic  and  polymer 
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characteristics  of  a  regular  composite  structure  is  to  introduce  dual  periods  rather  stricdy 
randcmi  periods.  This  would  be  simpler  to  fabricate  than  using  a  complete  distribution  of 
dimensions.  Figure  10  compares  the  attenuation  characteristics  of  a  regular  10%  structure, 
composed  of  0.2  mm  ceramic  strips,  and  a  structure  with  0. 1  and  0.3  mm  thick  ceramic 
strips.  The  pass-  band  still  exists  when  the  dual  period  geometry  is  employed  but  has  been 
fairly  effectively  broken.  Since  a  real  device  possesses  significant  losses,  which  are  not 
included  in  the  present  model,  one  would  expect  that  the  distortion  of  the  passband  of  a 
practkal  device  to  be  more  extreme. 

It  is  possible  to  achieve  a  similar  effect  by  filling  alternate  sawn  slots  with  polymers 
with  quite  different  elastic  properties.  In  practice,  every  second  slot  would  be  formed  and 
the  structure  filled  with  one  polymer.  In  a  second  operation,  the  second  sequence  of  slots 
would  be  formed,  and  these  filled  with  the  other  polymer.  Hence,  no  extra  sawing  is 
required  but  an  extra  polymer  cure  interval  is  added  to  the  manufacturing  process.  Figure 
1 1  illustrates  the  disruption  of  the  pass  band  achieved  in  10%  composite  with  0.2  mm 
ceramic  slots,  where  alternating  polymer  strips  with  wave  velocities  different  by  a  factor  of 
two  are  utilized. 

On  the  basis  on  the  limited  theoretical  work  presented  here,  it  is  evident  that  the 
amount  of  attenuation  encountered  is  a  function  of  volume  fraction,  variance  characteristics 
and  material  characteristics.  For  a  given  variance  fraction,  a  greater  attenuation  effect  is 
achieved  if  it  is  applied  to  the  phase  in  which  the  wave  spends  most  time  (i.e.  phase  for 
which  dimensicMi/velocity  is  greatest). 

(b)  Results  Observed  in  Practical  Transducers 

At  present,  theoretical  analysis  of  distributed  period  composites  is  limited  to  one 
dimensional  examples.  However,  experimental  tests  of  the  effect  of  employing  distributed 
periods  in  two  dimensicms  has  been  performed.  Some  one  dimensional  samples  were 
manuf^tured  but  a  more  visible  effect  was  observed  in  the  two  dimensional  samples. 
Furtheimore,  one  dimensitmal  compcsites  are  not  encountered  in  practice  to  any  noticeable 
extoit.  Hgure  12  illustrates  a  50%  ceramic  volume  fraction  composite.  Both  devices  have 
mean  cmmic  pillar  cross-sectional  dimensions  of  0.7  rrun,  but  the  distributed  dimension 
device  has  a  25%  variance  applied  to  the  ceramic  dimension,  across  the  structure  in  both  X 
and  Y  directions.  Figure  13  illustrates  tite  measured  inqjedance  magnitude  spectra  of  these 
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Figure  10 

Lateral  transmission  through  a  400  layer  structure 

Single  period  (0.2  ceramic,  1.8  polymer)  - 

Double  period  (0.1/0.3  mm  ceramic,  1.8  mm  polymer) - 


kHz 


Figure  11 

Lateral  transmission  through  a  400  layer  structure 
Regular  structure  (0.2  mm  ceramic,  1.8  mm  RE2039  polymer) 
Multi-poljrmer  (0.2  mm  ceramic,  1.8  mm  alternating  RE2039  polymer 
and  a  polymer  with  wave  velocity  half  that  of  RE2039) 
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rigure  13 

impedance  magnitude  of  50%  composite  transduc 
Single  period - - 

25%  ceramic  variance 


transducers.  A  reduction  in  spurious  resonant  activity  is  clearly  evident.  A  second  pair  of 
transducers,  with  a  10%  volume  fraction,  were  also  tested.  In  this  case  the  transducer  with 
a  25%  variance  in  the  polymer  dimension  also  demonstrates  reduced  spurious  resonant 
activity.  This  is  illustrated  in  Fig.  14. 

Figure  15  illustrates  results  obtained  for  25%  ceramic  volume  fraction  composites 
with  single  and  double-periodic  structures.  The  single  period  structiu-e  has  ceramic  lateral 
dimensions  of  1.0  mm,  and  the  double  period  structure  dimensions  of  0.4  and  1.2  nun 
alternately.  Significant  suppression  of  a  strong  spurious  mode  in  the  vicinity  of  the 
thickness  mode  is  evident 

It  has  been  theoretically  and  experimentally  demonstrated  that  distributed  period  and 
double  period  dicing  of  1:3  composites  is  effective  in  suppressing  spurious  (bandedge) 
resonances.  This  eliminates  the  need  for  very  frne  dicing  in  order  to  shift  the  modes  to 
higher  frequencies.  However,  overall  transducer  performance  must  still  be  optimized  by 
judicious  choice  of  ceramic  volume  fraction  and  mean  pillar  width-to-height  ratio. 

3.  0:3  COMPOSITE  TRANSDUCERS 

0:3  Composite  transducers,  comprising  microscopic  piezoceramic  particles 
dispersed  in  a  polymer  matrix,  have  proven  particularly  difficult  to  nxxlel  in  a  quantifrable 
manner.  As  discussed  in  the  interim  report  [1],  most  researchers  have  employed  the 
“cubes”  model,  or  a  “modified  cubes”  model  [4].  These  are  clearly  inadequate  since  the 
geometries  on  which  they  are  based  are  totally  unrelated  to  those  encountered  in  practice. 
Consequently  this  model  has  been  largely  disregarded  in  the  current  work. 

A  review  of  experimental  observations  regarding  the  dependence  of  0:3  transducer 
performance  on  manufacturing  details  has  been  made.  It  has  become  plainly  evident  that  the 
best  results  are  obtained  using  a  tight  packing  of  globular  grains  of  ceramic.  This  geometry 
implies  that  a  ceramic  volume  fraction  of  60  -  65%  be  used.  Devices  with  volume  fractions 
below  this  range  do  not  work  satisfactorily  because  of  the  presence  of  insulating  layers  of 
polymer  between  the  panicles.  It  is  physically  impossible  to  construct  a  0:3  composite, 
with  a  randcHn  packing  of  particles,  which  has  a  volume  fraction  greater  than  ^proximately 
65%.  Hence,  it  is  only  relevant  to  consider  this  narrow  practical  volume  range.  However, 
during  the  transition  frrom  a  suspension  of  particles  in  the  polymer,  to  a  structure 
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Figure  14 

Measured  impedance  magnitude  of  10%  composite  transducer 

Single  period  - 

25%  polymer  variance 
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Figure  15 

Measured  impedance  magnitude  of  25%  composite  transducer 

Single  period  - 

Double  period 


29 


comprising  an  interlocking  matrix  of  touching  particles,  the  net  properties  switch  from 
being  dominated  by  the  polymer's  properties  to  being  a  complex  function  involving 
principally  the  properties  of  the  ceramic  particles. 


3.1  New  Theory  for  the  Effective  Medium  Properties  of  Q:3_C! 


tes 


The  effective  medium  pr(^)erties  of  a  packing  of  spheres  is  a  function  of  the  number 
of  points  of  contact  and  the  conductance  of  each  contact.  The  area  of  particle  contact  is  a 
function  of  applied  pressure  and  the  elastic  properties  of  the  particles.  Equations  relating 
these  factors  have  been  combined  to  yield  one  equation  expressing  effective  medium 
pr(q)erties  (thermal  conductivity,  electrical  conductivity,  permitivity  etc.)  as  a  function  of  all 
of  these  factors. 

Nomenclature; 

a  Radius  of  particle 

r  Radius  of  circular  area  of  contact  between  adjacent  particles 
p  Number  of  points  of  contact  per  particle 

k  Electrical  conductivity,  or  thermal  conductivity,  or  dielectric  permittivity 
k*  Effective  ctniductivity  of  the  packing 

km  Conductivity  of  the  matrix 

kp  Omductivity  of  the  particles 

H '  Nondimensional  conductance  of  each  interface 

H  Dimensional  conductance  of  each  interface  (units  'S'  in  the  case  of  electrical 
conduction) 

V  Potential  across  interface 

E  Young's  modulus  of  elasticity  of  each  particle 

F  Force  pressing  adjacent  particles  through  the  axis  which  joins  their  centers 
Ph  Hydrostatic  (omnidirectional)  pressure  upon  packing 
0  Particle  volume  as  a  fractitm  of  total  packing  volume 
a  Poisson’s  ratio  of  the  particles 


In  the  case  of  a  piezoelectric  ceramic  powder/  polymer  cranposite,  the  ceramic  has  a 
far  higher  dielectric  permittivity,  electrical  conductivity  and  bulk  elastic  modulus  than  the 
polymer  phase.  Therefore,  when  any  of  the  simple  models  for  the  net  effective  properties 
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of  composite  materials  are  used  [5,6,7],  a  wide  gap  exists  between  the  upper  and  lower 
theoretical  bounds.  Experimental  evidence  suggests  that  the  actual  properties  do  not  fit  with 
any  of  the  limiting  bounding  cases.  Previous  modeling  of  this  type  of  transducer  had 
typically  involved  using  a  “cubes”  or  “modified  cubes”  model  [4].  This  involves  a  simple 
elastic  analysis  of  a  ceramic  cube,  or  cuboid,  surrounded  by  a  thin  film  of  polymer.  The 
aspect  ratio  of  the  cuboid  was  modified  until  a  match  between  theory  and  experiment  was 
achieved.  The  requirement  for  an  improved  strategy  is  readily  apparent.  The  approach 
presented  here  is  far  more  realistic.  It  is  based  on  a  quantitative  analysis  of  the  particle- 
particle  interactions  in  a  randomly  oriented  packing  of  isotropic  spherical  particles  of 
uniform  size.  Although  practical  piezoelectric  composites  are  not  composed  of  particles 
possessing  such  a  regular  geometry,  the  analysis  is  of  value  since  it  highlights  important 
factors  concerning  particle  properties,  the  number  of  points  of  contact  and  the  containing 
pressure.  Furthermore,  a  comprehensive  analysis  of  various  transducer  designs  indicates 
that  a  close  packing  of  approximately  spherical  particles,  with  a  narrow  size  distribution,  is 
a  highly  desirable  design  objective  [8,9,10].  Although  the  methodology  outlined  below  is 
capable  of  quantitative  analysis  when  all  the  constituent  parameters  are  known,  it  is 
probably  of  more  immediate  value  as  a  means  of  understanding  the  qualitative  relationship 
between  the  effective  medium  properties  and  the  various  factors  which  influence  these 
properties. 

It  is  known  that  the  physical  properties  of  electrical  conductivity,  thermal 
conductivity  and  dielectric  permittivity  are  analogous  [11].  Additionally,  the  bulk  elastic 
modulus  is  loosely  related.  It  is  not  related  exactly,  since  in  the  elastic  case  an  applied 
“potential”  (stress)  result  in  a  'flux'  (strain)  in  perpendicular  directions.  This  phenomenon 
does  not  appear  in  the  conductivity  and  permittivity  of  the  assumed  isotropic  media. 
Piezoelectric  effective  medium  theory  must,  of  ccHirse,  assume  anisotropic  particles. 

Batchelor  and  O’Brien  [12]  have  considered  the  problem  of  conduction  in  a  packed 
granular  structure  in  detail.  Their  technique  was  developed  for  isotropic  spherical  particles 
of  uniform  size.  Additionally,  it  was  necessary  to  assume  that  the  dielectric  permittivity  or 
conductivity  of  the  particles  was  far  greater  than  that  of  the  matrix  in  which  they  are 
embedded.  In  the  present  case  the  permittivity  of  the  polymer  is  approximately  fifty  times 
greater  than  that  of  the  polymer  and  the  conductivity  is  one  thousand  times  greater. 

Batchelor  obtained  an  expression  for  thermal  conductivity  as  a  function  of  particle 
volume  fraction,  the  number  of  points  of  contact  per  particle  and  the  “conductivity”  (non- 
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dimensional)  per  contact.  This  was  obtained  by  summing  the  effects  of  all  the  contacts  in  a 
specific  volume.  Exact  solutions  were  presented  for  geometrically  ordered  packings. 
However,  for  statistically  random  isotropic  geometric  structure  with  contact  points 
uniformly  distributed,  die  following  expression  was  obtained. 

^  =  (15) 

km  2 

Batchelor  asserted  that  if  the  flux  through  the  contact  area  between  touching 
particles  was  very  large  compared  with  that  through  the  adjacent  matrix,  then  the 
distributimi  of  potential  inside  the  particle  is  similar  to  that  found  in  the  velocity  potential  in 
irrotadonal  flow  of  incompressible  fluid  through  a  circular  hole  in  a  plate.  The  total  flux  in 
this  case  is  known  [13]. 

H  =  2  kp  V  r  (16) 

This  value  may  be  substituted  for  H'  in  Eq.  (15).  The  “non-dimensional”  flux  term  used  in 
Eq.  (15)  may  be  converted  into  a  real  “dimensional”  quantity  using  the  following 
transformation  presented  by  Batchelor  [12]. 

H  =  p  kp  V  a  H  (17) 

At  this  stage,  if  Eqs.  (15)  to  (17)  are  combined,  then  the  only  unknown  parameters 
are  the  radius  of  the  contact  surface  and  the  number  of  point  of  contact  per  particle.  The 
radius  of  the  contact  circle  is  a  function  of  the  force  pressing  adjacent  particles  together  and 
thdr  elastic  properties.  The  radius  of  the  contact  circle  may  be  evaluated  using  Hertz  contact 
theory  [14]. 


r  =  f1/3(  3  n-q2)  a  )l/3  (18) 

E 

Fortunately,  die  compressing  force  between  adjacent  particles  may  be  expressed  in 
terms  of  the  hydrostatic  (omnidirectional)  pressure  experienced  by  the  packing.  This  is 
obtained  by  dividing  the  total  force  experienced  by  each  particle  due  to  hydrostatic  pressure 
by  the  number  of  contact  points  and  modifying  the  result  to  take  account  of  finite  particle 
vdume  fraction  [15]. 
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(19) 


F  =  Ph4jLa2<> 
p 

Qxisequently,  Eqs.  IS  to  19  may  be  combined  to  yield  the  following  expression 
for  die  effective  medium  pix^ierty. 

k*  =  kp  (  (iEhUja^Dl/S  (20) 

7C  E 

From  the  analogy  noted  in  [1 1],  this  general  expression  may  be  used  to  calculate  the 
effective  electrical  conductivity,  thermal  conductivity  or  dielectric  permittivity  of  a  random 
packing  as  a  function  of  particle  volume  fraction,  the  number  of  points  of  contact  per 
particle,  the  applied  hydrostatic  pressure,  elastic  properties  of  the  particles  and  the  relevant 
“conductivity”  of  the  particle  phase.  However,  major  restrictions  still  apply  to  feasible 
values  for  <|>  and  p. 

In  mdered  packings,  the  particle  volume  fraction  may  vary  from  0.524  for  a  simple 
cubic  array  to  0.740  for  a  face  centered  cubic  array.  Also,  the  number  of  points  of  contact 
per  particle  varies  from  6  for  the  simple  cubic  case  to  12  for  the  face  centered  cubic  case. 
Although  the  range  of  possible  volume  fractions,  and  number  of  points  of  contact  appears 
wide,  the  range  obtained  for  packings  allowed  to  form  randomly  is  surprisingly  small. 
W.  O.  Smith  et  al  [16]  obtained  a  mean  volume  fraction  of  0.64.  He  used  a  packing  of  lead 
shot  which  was  soaked  in  acid.  When  the  acid  had  evaporated,  small  rings  of  lead  acetate 
(Ml  the  surfaces  of  the  spheres  revealed  the  l(x:ations  of  points  of  contact.  Upon  breaking  up 
the  packing  and  counting  these  rings,  the  authors  observed  that  the  mean  number  of 
contacts  per  sphere  was  9.  W.  O.  Smith  also  proposed  that  the  number  of  contacts  was  a 
linear  functi<Mi  of  v(dume  fracticni  between  the  theoretical  limits — i.e.,  6  contacts  per  sphere 
at  a  volume  fracticm  of  0.52,  to  12  points  of  contact  at  a  volume  fraction  of  0.74.  G.  D. 
Stxm  [17]  observed  that  the  volume  fraction  of  a  packing  of  uniform  spheres  varied  from 
0.60  to  0.64  depending  on  the  extent  of  agitation.  Bernal  and  Mason  [18]  observed  a  mean 
volume  fracticMi  of  0.62  using  a  technique  similar  to  that  of  W.  O.  Smith  et  al  [16],  but 
using  black  paint  instead  of  acid.  The  mean  number  of  points  of  contact  per  sphere  was 
found  to  be  approximately  six.  The  discrepancy  between  this  result  and  that  of  W.  O. 
Smith  is  periiaps  due  to  slightly  differing  definitions  of  what  “touching”  means  (i.e. 
physical  ccMitact  or  just  close  proximity).  Since  the  present  theory  assumes  that  physical 
pcmt  ccMitact  exists  befcMe  the  af^lication  of  hydrostatic  pressure,  the  lower  figure  for  the 
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number  of  contacts  will  be  chosen  here.  It  may  be  observed  that  while  particle  size  was  a 
factOT  in  most  of  die  preceding  equations,  it  is  not  in  Eq.  (20). 

It  is  reasonable  to  assume  that  the  actual  particle  volume  fraction  is  closely  related  to 
applied  hydrostatic  pressure.  Marion  et  al  [19]  have  performed  an  experimental 
investigatkm  of  the  relationship  between  particle  volume  fraction  and  applied  hydrostatic 
pressure  in  a  packing  of  silver  coated  glass  beads.  This  investigation  produced  the  curve 
illustrated  in  Fig.  16.  It  is  clear  that  for  volume  fractions  in  the  range  0.63  to  0.65, 
hydrostatic  pressure  and  particle  volume  fraction  are  linearly  related.  The  transition  from 
loose  suspension  to  tightly  bound  packing  is  remarkably  clear.  The  volume  fraction  at 
which  this  transition  occurs,  0.63,  matches  that  established  by  previous  workers  [16,17]. 
This  relatxMiship,  like  Eq.  (20),  is  only  valid  for  a  very  small  range  of  volume  fractions  just 
above  the  critical  packing  density  where  the  system  forms  a  packing  rather  than  a 
suspension.  Assuming  this  relationship  between  particle  volume  fraction  and  applied 
pressure,  the  volume  fraction,  in  Eq.  (20)  may  be  replaced  by,  cPh,  where  c  is  a 
constant  of  pnyortionality.  Hence,  the  following  pn^xntionality  expression  is  obtained. 

k*  »  c'Ph  (21) 


where  c'  is  a  constant  of  proportionality. 

The  effective  conductivity  is  a  linear  function  of  applied  hydrostatic  pressure. 
Marion  also  measured  the  electrical  conductivity  of  a  packing  of  silver  coated  glass  beads  as 
a  function  of  particle  volume  fraction  and  applied  hydrostatic  pressure.  The  result, 
illustrated  in  Fig.  17,  indicates  that  conductivity  is,  to  a  close  approximation,  a  linear 
function  of  applied  pressure.  It  should  be  noted  that  this  result  was  obtained  using 
insulating  beads  with  conducting  coatings.  Unfortunately,  it  was  impossible  to  measure  the 
conductivity  of  solid  conducting  beads  since  their  conductivities  are  extremely  high. 
Nevertheless,  the  conductivity  measured  here  is  dominated  by  the  particle-particle  contact 
mechanism  and  hence  the  result  obtained  still  serves  to  validate  the  theory.  Another 
convincing  example  of  the  validity  of  the  above  expression  was  the  observation  by  Edison 
[20]  that  the  electrical  conductivity  of  a  packing  of  carbon  spheres  is  a  function  of  applied 
jnessure.  This  characteristic  is  used  in  carbon  particle  telephone  microphones. 

The  accuracy  of  the  model  has  been  tested  in  a  limited  sense  by  considering 
publislwd  results  fra*  PbTiOs  particle  composites  [8].  In  this  case,  the  ceramic  volume 
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Figure  17 

Dependence  of  measured  electrical  conductivity  on  applied  hydrostatic  pressure  fOT 
a  packing  of  silver  coated  glass  beads.  Repn^uced  with  permission  of  A.  Nur. 
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firaction  was  approximately  60%  and  the  stmcture  was  formed  under  a  pressure  of  40  MPa. 
The  Young's  modulus  of  PbTiOs  is  150  x  10^  Nm‘2  and  its  Poisson's  ratio  is  0.20.  The 
relative  permittivity  of  the  pure  ceramic  is  approximately  230.  The  Voigt  model  predicts  a 
cmnposite  permittivity  of  140  while  the  Reuss  model  predicts  a  permittivity  of  12.  A 
relative  permittivity  of  23  is  predicted  using  Eq.  (6).  Although  the  observed  permittivity 
(45)  is  significantly  higher,  this  new  model  has  produced  a  closer  estimate  than  either  the 
Voigt  or  Reuss  models.  Furthermore,  no  arbitrary  “shape  factor”  has  been  employed,  as  is 
generally  used  in  the  conventional  “modified  cubes”  model  [7].  Ideally,  the  accuracy  of  the 
new  model  should  be  improved  by  attaining  a  better  understanding  of  the  number  of  contact 
points,  and  their  areas,  for  cases  where  particles  are  not  perfectly  spherical  and  do  have  a 
finite  size  distribution. 

An  expression  for  determining  the  conductivity  of  a  close  random  packing  of  high 
conductivity  spheres  has  been  developed.  It  takes  account  of  all  the  key  factors.  These 
include  particle  conductivity,  particle  elastic  properties,  applied  pressure,  particle  volume 
fraction  and  the  number  of  points  of  contact  per  particle.  The  expression  is  valid  only  for  a 
very  limited  range  of  particle  volume  fractions.  However,  this  range  encompasses  those 
which  are  of  most  practical  importance.  Using  an  observed  linear  relationship  between 
particle  volume  fraction  and  applied  pressure,  the  expression  correctly  predicts  the  widely 
observed  linear  relationship  between  effective  conductivity  and  applied  pressure.  The 
exixession  gives  a  limited  insight  into  the  parameters  critical  to  the  successful  manufacture 
of  piezoelectric  ceramic  powder  composites.  However,  it  is  known  that  the  physical 
mechartisms  governing  the  permittivity  in  ferroelectric  powders  are  rather  complex  [21]. 
For  example,  the  permittivity  near  the  surface  of  a  ferroelectric  particle  is  less  than  that  in 
the  interior.  Additionally,  the  permittivity  and  piezoelectric  coefficients  of  practical 
ferroelectric  powders  are  a  function  of  particle  size.  The  technique  employed  here  has 
assumed  uniform  properties  independent  of  size.  Incorporation  of  these  phenomena,  and 
development  of  effective  medium  models  for  elasticity  and  piezoelectricity,  will  require 
furtiier  exploitation  of  the  principles  presented  here. 

3.2  Observations  Made  Regarding  Factors  Crucial  to  0:3  Composite  Manufacture 

The  technical  literature  relating  to  0:3  composite  transducer  design  is  now 
extensive.  A  shmt  review,  collecting  generally  observed  factors  relevant  to  successful 
manufactioe,  is  presented. 
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The  observaticMi  that  a  ceramic  volume  fracticm  in  the  range  60%  to  70%  is  desirable 
is  general  [8,9^2,23].  Additionally,  the  best  results  have  been  obtained  using  particles 
which  are  rounded  in  shape  and  which  possess  a  narrow  distribution  of  size  [10.24,25] 
Two  techniques,  which  give  similar  net  results,  are  available.  The  co-precipitation 
technique  developed  by  Lee  et  al  [10],  involves  mixing  Pb(N03)2,  TiCw,  H2O2  and 
NH4OH  into  a  solution.  The  precipitate  is  dried  at  100°  C  and  calcined  at  800-900°  C  to 
produce  globular  particles  with  a  size  of  around  4pm.  The  more  conuTX>n  route  to  obtaining 
desirable  particle  characteristics  involves  quenching  the  ceramic  immediately  after 
calcination.  The  ceramic  breaks  down  easily  to  virtually  single  domain  crystallites  due  to 
the  high  strain  ratio  found  in  PbTiOs  [25].  Composites  manufactured  using  powder 
obtained  by  milling  solid  ceramic  have,  without  exception,  displayed  inferior  piezoelectric 
properties  [22]. 

Figure  18  illustrates  the  particle  size  distribution  obtained  using  the  quenching 
technique  [25].  Fuithermore,  the  size  of  the  particles  may  be  controlled  by  altering  the  PbO 
content  of  the  powder  during  calcination  [25].  Giniewicz  [24]  also  states  that  particle  size 
may  be  controlled  by  varying  firing  temperaturc/time,  or  by  drying  the  Pb'n03  with  Bi  and 
Fe.  NTK  further  improved  the  quality  of  their  powder  by  etching  it  in  an  acid  (eg  HCl). 
This  has  the  effect  of  dissolving  very  small  particles  (outside  the  desired  size  range), 
rounding  any  sharp  particle  comers  and  leaving  the  surfaces  well  prepared  for  wetting  to 
the  polymer  [25]. 

There  is  general  agreement  that  PbliOs  should  be  used  as  the  ceramic  phase.  This 
is  because  of  its  low  permittivity  and  high  %yi  compared  to,  for  example,  PZT.  However, 
several  different  polymers  have  been  used  successfully.  NTK  use  “Fluororubber,”  as 
described  in  detail  in  their  patent  [25].  Giniewicz  [24],  who  has  produced  composites 
possessing  dhgh  coefficients  comparable  to  those  of  NTK,  recommends  the  use  of  a  soft 
polymer  cm  the  basis  that  it  might  better  accommodate  realignment  of  the  domains  during 
poling.  Low  viscosity  is  obviously  desirable  and  Eccogel  1365-25  [26]  has  been  used  with 
success.  Giniewicz  recommends  a  thorough  technique  for  mixing  the  phases  [24].  The 
filler  is  passed  through  a  sieve  directly  into  the  prepared  polymer  to  prevent  unnecessary 
aggkMnmtion.  The  mixture  is  then  shear  blended  using  a  broad  flexible  knife  to  maximize 
die  degree  of  wetting  between  the  coamic  and  polymer. 
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Figure  18 

Particle  size  distribution  obtained  by  NTK  (“percentage  particles  smaller  than”) 

(a)  PbTiOa  ceramic  quenched  directly  fnan  oven. 

(b)  PbTiOa  ceramic  flowed  to  cool,  then  milled. 

(c)  P27r  ceramic  quenched  direaly  from  oven. 

(d)  PZT  ceramic  ^owed  to  cool,  then  milled. 
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The  pressure  used  during  formation  of  the  composite  prior  to  curing  is  also  critical. 
Westinghouse  [27]  typically  used  low  pressures  of  around  6  MPa.  It  is  believed  that  their 
use  of  low  pressure  during  formation,  combined  with  the  use  of  less  than  perfect  ceramic 
particles,  contributed  to  their  major  poling  problems.  NTK  [25]  claim  to  use  a  pressure  of 
approximately  14  MPa.  Amongst  the  numerous  results  obtained  by  Pennsylvania  State 
University,  the  best  dhgh  results  were  obtained  when  the  highest  reported  pressure  was 
used.  Lee  et  al  [10]  report  obtaining  a  dhgh  coefficient  of  4170  x  10'15  m^N-l  when  a 
pressure  of  65  MPa  was  used  to  compress  globular  4  pm  particles  obtained  by  the  co¬ 
precipitation  method. 

After  curing  under  pressure  according  to  the  temperature/time  schedule 
recommended  for  the  polymer,  the  composite  is  ready  for  electroding.  The  surfaces  are 
generally  lightly  lapped  to  expose  a  bare  ceramic  surface.  Generally,  silver  paint  has  been 
used  [8,10]  but  NTK  have  successfully  used  a  polymer  loaded  with  carbon  black,  as 
described  in  detail  in  their  patent  [25].  A  carbon  loaded  polymer  (similar  polymer  to  that 
within  the  composite)  may  have  the  advantages  of  greater  robustness  and  temperature/ 
expansion  compatibility  with  the  composite  material. 

Generally,  poling  has  been  performed  under  the  following  conditions.  A  DC  field 
of  75- '  JJ  kVcm'^  is  applied  for  approximately  30  minutes  while  the  sample  is  immersed  in 
oil  w  75°C-1000C  [8,10].  The  field  should  probably  be  ramped  up  and  down  over  a  short 
time  interval.  Since  the  transducer  is  essentially  a  capacitor,  rapid  fluctuations  in  voltage 
will  result  in  current  surges  which  do  not  help  to  align  the  domains  but  may  result  in 
damage  in  the  form  of  electrical  breakdown  and/or  thermal  runaway.  Corona  poling  has 
also  been  tried  [28],  but  when  directly  compared  with  results  obtained  by  conventional 
poling,  the  difference  in  d33  and  Ef  is  scarcely  significant.  However,  in  the  case  of  corona 
poling,  if  breakdown  occurs  at  one  location,  then  poling  of  the  rest  of  the  sample  may 
continue.  In  conventional  poling,  any  breakdown  will  immediately  short  the  entire  high 
voltage  electrode  to  ground  and  a  destructive  current  surge  will  occur. 

Various  techniques  have  been  attempted  for  improving  poling  efficiency.  Sa-Gong 
et  al  [29]  considered  the  effect  of  adding  carbon,  germanium  and  silicon  fillers  to  the 
polymer.  Since  the  conductivity  of  the  ceramic  phase  is  several  orders  of  magnitude  higher 
than  that  of  the  polymer  phase,  there  is  a  danger  that  most  of  the  applied  poling  field  will  be 
across  the  polymer  phase.  Typically  the  poling  field  used  for  a  0:3  composite  is  five  times 
higher  th  n  that  which  would  be  used  across  pure  ceramic.  Therefore,  if  there  are 
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inconsistencies  (eg  agglomerations)  in  the  composite  structure,  then  localized  fields  which 
occur  may  be  sufficient  to  cause  breakdown.  The  addition  of  semiconducting  fillers  did 
ease  poling  but  resulted  in  higher  dielectric  loss  in  the  resulting  material.  However,  the 
best  results  obtained  were  inferior  to  those  obtained  by  Lee  et  al  [10]  who  used  no  fillers. 
Giniewicz  considered  improving  poling  performance  by  doping  the  PbTiOs  with  Mn  [24]. 
This  has  the  effect  of  lowering  the  conductivity  of  the  ceramic,  and  hence  enhancing  the 
field  across  the  particles  during  poling.  An  improvement  of  40%  in  dhgh  figure  of  merit 
was  obtained.  However,  even  this  result  is  lower  than  that  obtained  by  Lee  et  al  [10].  In 
fact  the  d33  obtained  by  Lee  is  virtually  identical  to  that  observed  in  single  phase  PbTi03. 
This  implies  that  almost  ideal  polarization  has  been  achieved.  It  is  concluded  that  Lee’s 
method  [10]  is  the  best.  Fortunately  it  is  simple  in  that  no  special  fillers  are  used  and  relies 
upon  very  high  pressure  formation  using  globular  particles  with  a  narrow  size  distribution. 

Another  process  worthy  of  mention  is  the  “fired  composite”  technique  developed  at 
Pennsylvania  State  University  [28,30,31],  PbTi03  powder  is  prepared  by  the  co¬ 
precipitation  technique  described  by  Lee  [10].  Therefore,  it  can  be  assumed  that  the 
particles  are  globular  and  of  narrow  size  distribution.  The  ceramic  powder  was  mixed  with 
PVA  binder  and  pressed  into  pellets  at  extremely  high  pressure  (140  MPa).  The  pellets 
were  then  placed  in  an  oven  and  fired  at  temperatures  of  up  to  900°C.  The  best  results,  in 
terms  of  highest  d33,  were  obtained  when  the  highest  temperature  was  used  [31].  In  one 
case  corona  poling  was  employed  and  found  to  produce  marginally  better  results  than 
conventional  poling  [28].  In  this  case  poling  was  performed  prior  to  polymer  impregnation. 
It  was  stated  that  corona  poling  was  ineffective  when  performed  on  impregnated 
composites.  It  is  hypothesized  that  a  limited  amount  of  breakdown  occurred  during  corona 
poling  due  to  the  presence  of  air  voids.  This  may  have  the  beneficial  effect  of  placing  more 
of  the  applied  field  across  the  individual  particles.  Presumably,  if  conventional  poling 
techniques  are  applied  to  unimpregnated  composites,  catastrophic  breakdown  will  occur 
due  to  the  rapid  establishment  of  a  low  resistance,  high  current  path  between  the  high 
potential  top  electrode  and  the  grounded  electrode.  Composite  transducers  manufactured 
using  the  “fired”  technique  generally  have  high  volume  fractions  ranging  from  65%  to  70% 
— implying  extremely  compact  particle-particle  contact.  Consequently,  they  could  be  poled 
at  slightly  lower  fields  (50  kVcm'*).  The  highest  dhgh  figure  of  merit  reported  using  the 
“fired”  technique  was  1750  x  10-15  m^N'i,  which  is  less  than  half  of  that  obtained  by  Lee 
[10].  However,  this  figure  does  not  relate  to  a  device  formed  at  the  highest  (optimum) 
temperature. 
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The  key  characteristics  relevant  to  successful  manufacture  have  now  been 
established.  This  section  will  conclude  with  a  brief  an^ysis  of  the  importance  of  particle 
size.  Relatively  little  quantitative  comparable  information  is  available.  However,  NTK 
PR306  and  NTK  PR305  have  been  compared  by  Banno  [32].  PR306  is  composed  of 
particles  approximately  1  |im  across  while  PR305  is  composed  of  10  4m  particles  [33]. 
The  large  particle  material  (PR305)  has  a  dhgh  figure  of  merit  of  5084  x  lO'iS  m^N-l 
whereas  the  small  particle  material  (PR306)has  coefficient  of  1 160  x  lO'lS  m^N'^.  The  low 
pressure  reception  sensitivity  of  the  larger  particle  material  is  7dB  higher  than  that  of  the 
other  material.  However,  while  the  reception  sensitivity  of  the  small  particle  material  is 
virtually  independent  of  hydrostatic  pressure,  the  gh  of  the  large  particle  material  falls  by 
approximately  25%  as  the  pressure  is  increased  to  10  MPa  [32].  During  the  same  pressure 
increase,  the  relative  permittivity  of  the  large  particle  material  rises  by  10%  whereas  the 
permittivity  of  the  small  particle  material  does  not.  It  has  been  proposed  that  that  the 
pressure  instability  of  the  large  particle  material  is  due  to  voids  which  collapse  under  high 
pressure  [34].  Since  it  is  generally  desirable  for  the  sensitivity  of  a  hydrophone  to  be 
pressure  independent  and  tolerant  of  multiple  pressure  cycling,  the  conclusion  is  drawn  that 
smaller  panicles  are  preferable.  Giniewicz  [8]  obtained  reasonably  pressure  independent 
results  using  5  pm  particles.  The  composites  were  pressed  at  40  MPa  but  allowed  to  cure  at 
atmospheric  pressure.  An  interesting  observation  may  be  made  regarding  the  parameter 
pressure  dependency  illustrated  in  Fig.  19.  As  soon  as  pressure  is  applied  a  small  jump  in 
gh  sensitivity  is  observed.  It  is  hypothesized  that  the  particles  are  being  pressed  into  better 
contact.  The  sensitivity  declines  very  slightly  with  further  pressure  increase.  Since  dh  is 
constant  over  this  range,  the  slight  decrease  in  gh  must  be  due  to  increasing  permittivity, 
which  might  also  be  related  to  particle  contact.  On  release  of  pressure,  gh  is  generally  lower 
than  on  the  pressure  increase  cycle.  This  may  be  due  to  breaking  particle-particle  bonds. 
This  data  would  suggest  that  high  pressure  should  be  used  during  formation  and  curing  so 
that  subsequent  pressure  cycling  is  negligible  in  comparison  to  that  endured  during 
manufacture. 

In  summary,  the  optimal  results  are  obtained  using  PbTi03  powder  formed  by 
quenching  or  co-precipitadon.  The  ceramic  particles  are  globular  in  shape  and  have  a 
narrow  size  distribution.  The  particles  should  be  etched  in  acid  to  clean  them,  remove  sharp 
comers  and  to  eliminate  small  particles.  The  particles  and  the  chosen  polymer  (Eccogel 
1365-25  is  adequate)  are  carefully  blended  and  formed  into  tablets  at  high  pressure  (65 
MPa)  and  cured.  ConventitMial  poling  in  oil  at  75oC,  in  a  field  of  75-100  kVcm**  for  30 
minutes  is  sufficient  to  obtain  charactoistics  cotiqjetitive  with  the  best  available  [10,32] 
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FigiiTC  19 

Hydrostatic  Piezoelectric  Coefficients,  Observed  by  Giniewicz,  as  a 

Function  of  Pressure 


4.  TRANSDUCERS  COMPRISING  MULTIPLE  ACTIVE  LAYERS 


Nomenclature 

F  Force  ( N ) 

h  Piezoelectric  constant  ( N/C  or  V/m ) 

Co  'Qamped'  capacitance  of  transducer  ( F ) 

ti  Wave  transit  time  across  layer  'i'  ( s ) 

Vs  Voltage  across  source  ( V ) 

Vi  Voltage  across  electrodes  of  transducer  layer  'i' 

Zs  Electrical  impedance  of  source  ( Q ) 

Zi  Mechanical  impedance  of  layer  'i*  ( kgs'  ^ :  pvA ;  A  =  area ) 

Zp  Mechanical  impedance  of  front  load 

Zb  Mechanical  impedance  of  back  load 

s  Laplace  operator  (Bar  symbol  denotes  Laplace  quantities). 


Thickness  mode,  backed  and  matched,  piezoelectric  transducers  perform  a  central  role 
in  medical  ultrasound,  non-destructive  examination  and  high  frequency  sonar  applications. 
However,  transducer  performance  is  frequently  a  factor  limiting  ultimate  system  performance. 
The  concept  of  using  multiple  active  layers  provides  the  potential  for  simultaneously  enhancing 
both  bandwidth  and  sensitivity.  It  offers  the  prospect  of  optimized  transducer  performance,  and 
takes  advantage  of  recent  hardware  developments  that  have,  in  many  ways,  resulted  in  an 
imbalance  between  transducer  and  processing  hardware  technologies. 

Conventionally,  the  thickness  mode  transducer  has  been  composed  of  a  single 
piezoelectric  layer  and  one,  or  more,  matching  layers  serving  as  a  passive  piezoelectric 
transfcnmer.  It  has  recently  been  proposed  that  the  matching  layer  should  also  be  piezoelectric 
and  that  voltage  functions  be  applied  to  it  in  order  to  further  enhance  the  sensitivity  and 
bandwidth  of  the  transducer  [35,36].  The  configuration  of  the  new  transducer  is  illustrated  in 
Fig.  20.  Of  particular  importance  is  the  ability  to  design  a  transducer  with  a  frnite  output  at  a 
frequency  where  normally  a  null  would  exist.  Previous  work  [35,36]  concentrated  on  the 
transmissimi  characteristics  of  multiple  layer  transducers  in  which  independent  stimulation  was 
^lied  to  each  layer.  These  techniques  are  briefly  reviewed.  This  section  concentrates  on  the 
design  of  multiple  layer  transducers  as  receivers.  The  performance  of  a  multiple  layer 
transducer  in  the  reception  mode  is  critically  dependent  on  optimal  design.  However,  the  same 
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Layer  1  Layer  2 


Layer  1  Layer  2 


Figure  20 

Configuration  of  the  new  double  active  layer  transducer.  The  particle  displacement 
fuiKtions  within  each  layer  are  defined. 
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benefits  of  enhanced  sensitivity  and  bandwidth  are  achievable.  Hence,  high  performance  pulse- 
echo  c^)eration  is  feasible.  Additionally,  the  stability  of  the  technique  with  respea  to  deviations 
in  actual  design  parameters  is  investigated. 

4.1  Transmission 

Two  configurations  are  considered.  In  the  “active  matching  layer”  scheme  the  input  u> 
the  principle  active  layer  (layer  1)  is  pre-defined  A  L^lace  domain  matrix  technique,  based  on 
that  described  by  Lewis  [37],  has  been  used.  For  the  case  of  a  two  layer  system,  where  the 
excitation  voltages  are  specified,  one  obtains  the  matrix  problem  illustrated  in  Fig.  21  [36]. 
Notice  that  the  known  quantities  (input  voltages)  are  on  the  right  hand  side  and  the  known 
quantities  (output  displacements)  are  on  the  left  hand  side.  In  the  “active  matching  layei^’  case, 
the  output  force  and  the  excitation  to  layer  1  is  defined.  The  voltage  required  across  layer  2,  in 
order  to  satisfy  the  physical  constraints  on  the  system,  must  be  established.  The  matrix  is 
altered  to  obtain  the  solution  to  this  problem  by  exchanging  the  positions  of  the  relevant 
quantities  in  the  matrix  (enclosed  in  boxes  in  Fig.  21).  It  should  also  be  noted  that  by  placing  a 
quantity  pertaining  to  the  incident  force  on  the  right  hand  side,  it  is  possible  to  obtain  the 
received  t^n  circuit  voltage  response  across  each  layer  [36]. 

In  the  “(^timized  excitation”  technique  there  is  no  need  to  pre-specify  the  excitation  of 
any  layer.  Only  the  required  output  force  is  specified.  Although  the  result  is  a  multiple  input, 
single  output  problem,  with  an  infinite  number  of  solutions,  an  optimized  solution  is 
established. 

Since  we  are  dealing  with  a  linear  system,  the  total  force  output,  in  the  Laplace  domain, 
may  be  exfnessed  as  follows 

F  =  Pi  Vi  +  P2V2  + .  PnVn  (22) 

where  “n”  is  the  number  of  layers,  and  the  pi's  are  the  transfer  functions  between  the  applied 
Vintages  and  the  output  force. 

The  quantity  pi  may  be  obtained  by  measuring  F  when  Vi  =  1.0  and  Vi  =  0.0,  i »  2  to 
n.  Pi ,  is  2  to  n,  may  be  evaluated  in  a  similar  manner.  If  one  considers  a  circuit  analogy 
model  for  the  transducer  system,  this  is  simply  a  use  of  the  principle  of  current  superposition. 
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Ai  and  B]  etc.  are  coefficients  of  panicle  displacement  as  described  in 
Fig.  1 


Figure  21 

The  governing  matrix  for  the  double  active  layer  transmitter  example.  Known 
quantities  (inputs)  are  on  the  right,  and  unknown  quantities  (outputs)  are  on  the 
left.  If  requin^,  the  matrix  elements  may  be  rearranged. 
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The  net  current  in  the  output  loop  is  the  sum  of  those  due  to  the  voltage  sources  considered 
independently. 

For  the  sake  of  sin^licity,  let  Vi  =  1.0  /jO®  ,  i  =  1  to  n.  Then,  Eq.  22  reduces  to 

F  =  3i+P2  + . Pn  (23) 

or 

IF!  46  =  Ipll  401  +  Ip2l  402  + . +  Ipn'  (24) 

At  this  point,  the  discussion  will  be  limited  to  the  two  layer  case.  The  phasor  quantities, 
Pi  and  p2>  represented  graphically  in  Fig.  22.  Force  output,  IFI40  is  maximized  if  02  = 
01.  In  this  case  plVi  and  P2V2  aiv  in  phase  and  form  a  straight  line  when  summed.  This  may 
be  achieved  by  adjusting  the  phase  angle  of  V2  to  ctxnpensate  for  the  phase  difference  between 
Pi  and  P2.  The  desired  result  is  obtained  by  setting  V2=Vi  4(01-02)-  Hence,  for  this 
configuration,  Eq.  24  becomes 

«FmI401=  (iPlI  401  XIV1I4  0®)+(iP2l  402)OViI  4(01  -  ©2))  (25) 

where  Fm  is  the  maximized  value  of  output,  F.  Equation  25  represents  the  optimized  'transfer 
function*  of  the  system. 

Suppose  that  the  desired  output  is  IFI40'.  This  is  obtained  by  multiplying  Eq.  25 
throughout  by  IFi4.0’  . 

IFMU-ei 

Equation  25  then  becomes 

IFI40’  =  (lpll40l)(lVil4O®)(lEli^  +  (iP2J402)(lVll4(0l-02))(!E!4fl: ) 

IFMl4ei  IFMl4ei 


V’l 


V’2 


(26) 


In  this  way,  the  excitation  functions,  V'l  and  V’2  are  obtained.  Their  time  domain 
counterparts,  Vi(t)  and  V'2(t),  are  obtained  using  an  inverse  FFT. 
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Figure  22 

Phasor  diagram  illustrating  Pi  and  P2>  and  the  effect  of  phasing  P2  by 
multiplying  it  by  1.0  4(6i  -  O2)  and  adding  to  Pi. 
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In  this  way,  the  excitation  functions,  V'l  and  V'2  are  obtained.  Their  time  domain 
counterparts,  V'i(t)  and  V'2(t),  are  obtained  using  an  inverse  FFT. 

The  nature  of  the  optimized  “transfer  function”  of  a  double  layer  transducer  is  a 
function  of  its  design.  Clearly,  a  wide  smooth  response  is  desirable.  Figure  23  illustrates  the 
optimized  “transfer  functitm"  of  Transducers  A  and  B.  (The  transducer  designs  are  described  in 
Appendix  II.)  The  responses  obtained  for  the  cases  where  the  phase  angle  between  the 
excitation  functions  of  layers  1  and  2  were  kept  equal  to  zero  are  shown  for  comparison.  In  the 
case  of  both  transducers,  a  significant  output  is  achievable  at  the  second  harmonic  frequency 
simply  by  using  the  optimized  phasing  technique.  It  is  readily  evident  that  a  smoother,  and 
more  desirable,  response  has  been  obtained  using  Transducer  B,  which  is  more  heavily 
acoustically  dampened.  If  one  considers  that  both  sensitivity  and  bandwidth  are  quantities 
which  arc  to  be  maximized,  then  one  can  use  the  integral  of  sensitivity  with  respect  to 
frequency  as  a  useful  tigure  of  merit.  When  the  integrals  of  the  curves  illustrated  in  Figiue  23 
were  calculated,  it  was  found  that  by  using  optimum  phasing,  as  opposed  to  zero  phasing,  the 
figures  of  merit  of  both  Transducers  A  and  B  were  increased  by  approximately  90%. 

The  transmission  characteristics  of  Transducer  A  were  measured  experimentally.  The 
response  was  measured  in  the  far  field  (range  *  200  mm)  using  a  sensitive  hydrophone  [38] 
possessing  a  hydrophone/  pre-amplifier  sensitivity  of  -222  dB  (0  dB  =  1  V/pPa).  For 
comparison  with  experiment,  it  was  necessary  to  include  the  effect  of  diffraction  on  the 
theoretical  pressure  response  emanating  from  the  transducer  face.  The  transducer’s  diffraction 
impulse  response  was  calculated  for  this  geometry  [39]  and  convolved  with  the  pressure 
response  at  the  transducer  face.  In  this  case,  at  this  range,  the  diffraction  effect  is 
approximately  equivalent  to  differentiating  with  respect  to  time.  Consequently,  high  frequency 
components  are  emphasized.  The  result,  for  both  optimal  and  zero  phasing,  illustrated  in  Figs. 
24  and  25  indicates  that  the  transducer  is  performing  as  predicted. 

The  phase  angles  (0i  -  02)  required  to  achieve  the  optimal  “transfer  functions”  are 
illustrated  in  Fig.  26.  In  both  cases,  the  characteristics  are  smooth  with  no  discontinuities.  The 
phase  characteristic  of  Transducer  B  is  less  wavy  and  may  be  more  closely  approximated  by  a 
straight  line.  This  is  signifrcant,  since  a  straight  line  phase  characteristic  corresponds  to  a 
'delay'  in  the  time  domain.  A  transducer  possessing  a  perfectly  straight  optimum  phase 
characteristic  would  be  perfectly  matched.  This  implies  a  lack  of  reverberation  within  the 
transducer  structure.  While  it  is  readily  possible  to  compensate  for  reverberation  in 
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Rgure  23 

Transmission  “transfer  functions”  of  Transducer  A  and  Transducer  B. 
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Figure  24 

Transmission  “transfer  function"  of  Transducer  A,  using  optimal  phasing,  taking 
account  of  diffraction  effects. 


Figure  25 

Transmission  “transfer  function”  of  Transducer  A,  using  zero  phasing,  taking 
account  of  difliaction  effects. 
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Figure  26 

Phase  characteristics,  ( 0i  -  02  ),  to  be  applied  to  Transducer  A  and  Transducer  B 
in  order  to  achieve  the  optimal  “transfer  function.” 
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transmission,  by  canceling  its  effects  using  the  stimulation  functions,  this  is  not  possible  in  the 
reception  mode.  The  desirability  of  a  well  matched  transducer  structure  is  evident 

The  “optimum  excitation”  technique  has  been  tested  experimentally  for  non  harmonic 
signals  using  Transducer  C.  A  1.5  MHz  single  cycle  sinusoidal  pulse,  with  peak  amplitude  10 
kPa,  was  desired.  Using  the  above  technique,  the  required  excitation  functions  for  layers  1  and 
2  were  determined.  These  are  illustrated  in  Fig.  27.  Notice  that  they  are  simple  in  form  and  are 
similar  except  that  the  stimulation  for  layer  2  is  delayed  by  approximately  0.15  ms  with  respect 
to  layer  1.  The  waveforms  were  sampled  at  20  MHz  and  programmed  into  a  HP  8175A 
Arbitrary  Function  Generator  (AFG).  The  sensitivity  of  the  hydrophone  was  excellent  and 
hence  the  transducer  could  be  driven  directly  from  the  AFG  output.  The  response  was 
measured  at  a  range  of  3  mm  so  that  the  plane  wave  output  could  be  obtained.  This  response, 
illustrated  in  Fig.  28,  indicates  a  good  match  with  the  desired  waveform.  The  “tail,”  that  may 
be  observed  at  the  end  of  the  pulse,  has  been  attributed  primarily  to  an  imperfect  bond  layer.  It 
is  confidently  expected  that  transducers  with  substantially  better  bond  layers  can  be  made  in  the 
future.  Alternatively,  if  the  characteristics  of  the  bond  layer  can  be  quantified  comprehensively, 
ctxnpensation  for  them  can  be  made  when  the  excitation  functions  are  being  calculated. 

4.2  Reception 

As  in  the  transmission  case,  two  techniques  are  available.  The  first  of  these,  the  “active 
matching  layer,”  is  the  corollary  of  the  corresponding  technique  described  for  the  transmission 
case.  In  principle,  if  there  is  an  incoming  pressure  wave,  a  signal  may  be  applied  to  layer  2  so 
that  the  signal  monitored  across  layer  1  fits  some  pre-defined  function.  However,  the  technique 
is  of  limited  practical  importance  since  the  precise  nature,  amplitude  and  time  of  the  input 
function  must  be  known  in  advance.  Generally,  this  requirement  is  not  fulfilled. 

The  “delay  and  sum”  technique  is  related  to  the  “optimal  excitation”  method  used  in 
transmission.  This  technique  has  also  been  described  recently  by  Chofflet  et  al  [40].  It  has  been 
found  that  in  a  perfectly  matched  transducer,  the  optimum  phase  characteristics,  describing  the 
excitation  of  layer  2  with  respect  to  layer  1,  was  closely  related  to  the  transit  delay  of  a  wave 
propagating  through  layer  2.  Since,  in  the  reception  case,  the  inputs  are  time  domain  signals,  it 
is  desirable  to  approximate  the  exaa  phase  requirement  with  a  simple  time  delay  which  should 
be  applied  to  the  response  of  layer  2  prior  to  summing  with  the  response  of  layer  1.  Using  this 
technique,  a  finite  sensitivity  can  be  achieved  continuously  from  below  the  fundamental 
rcs(Miant  frequency  to  above  the  third  harmonic.  However,  the  “delay  and  sum”  technique  does 
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Figure  27 

Excitation  functions  to  be  applied  to  Transducer  C  to  achieve  a  1.5  MHz,  single 
cycle  pressure  profile  with  amplitude  10  kPa 


Figure  28 

Experimentally  observed  waveform  when  the  functions  illustrated  in  Fig.  27  were 
applied  to  Transducer  C. 
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not  eliminate  the  effects  of  reverberation  within  the  transducer.  Reverberation  is  only  a 
significant  problem  at  higher  frequencies  where  the  wavelength  is  short  relative  to  the  layer 
dimensions.  The  solution  to  the  problem  of  reverberation  is  to  use  a  perfectly  matched  backing 
to  eliminate  reflection  at  the  layer/  backing  interface. 

Figure  29  illustrates  the  receiver  voltage  output  of  Transducer  D  using  the  “delay  and 
sum”  technique.  Responses  are  shown  for  a  two  layer  system  for  the  case  where  a  0.13  |is 
delay  was  applied  to  the  response  of  layer  2  and  for  the  case  where  no  delay  was  applied.  The 
substantial  improvement  in  usable  bandwidth  is  readily  evident.  Alternatively,  if  we  consider 
the  response  of  a  single  layer  transducer,  with  thickness  half  that  of  the  two  layer  system,  one 
obtains  a  similar  bandwidth  to  the  “delay  and  sum”  result,  but  with  only  half  the  sensitivity. 

Figure  30  illustrates  the  response  of  Transducer  D  obtained  for  three  Gaussian  pulses 
with  center  frequencies  of  1.8,  3.6  and  5.4  MHz  and  peak  amplitudes  of  10  kPa.  These 
correspond  approximately  to  the  frequencies  of  the  fundamental,  second  and  third  harmonics. 
The  responses  have  been  displayed  so  that  they  are  not  superimposed.  Undistorted  pulses  have 
been  obtained  in  each  case.  The  response  obtained  when  using  Transducer  C  to  monitor  the  5.4 
MHz  signal  is  shown  for  comparison.  In  this  case  some  distortion  is  evident  due  to 
reverberation  in  the  active  layers.  Although  it  is  not  necessary  to  have  matched  backing  in  order 
to  obtain  finite  output  by  using  the  “delay  and  sum”  technique,  it  is  necessary  if  the  effects  of 
reverberation  on  high  frequency  signals  are  to  be  totally  eliminated.  Unfonunately,  practical 
multiple  layer  transducers  have  not  been  tested  in  the  reception  mode.  Since  the  transducers 
operate  largely  according  to  theoretical  prediction  in  transmission,  it  is  reasonable  to  assume 
that  their  practical  reception  performance  also  matches  theoretical  results. 

4.3  PulscrEcho  AppUcanons 

It  is  clearly  possible  to  combine  a  transmission  system  using  the  “optimized  excitation” 
technique  with  a  reception  system  using  the  “delay  and  sum”  technique.  As  discussed  above,  it 
is  preferable  to  use  a  matched  backing  in  order  to  optimize  the  reception  performance.  If 
piezoelectric  ceramic/  epoxy  composites  or  piezoelectric  polymers  are  employed,  this  constraint 
is  not  unduly  burdensome.  The  electronic  hardware  which  would  be  required  to  achieve  the 
required  functions  for  the  transmission  mode  would  not  need  to  be  particularly  elaborate. 
Synthetic  waveforms  may  be  stored  digitally  in  memory.  Sequential  digital  values  may  be  read 
out  and  converted  to  analog  signals  prior  to  amplification.  The  “delay  and  sum”  function 
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Rgure  29 

Theoretical  reception  magnitude  characteristic  of  Transducer  D. 


Figure  30 

Theoretical  reception  response  of  Transducer  D,  using  the  “delay  and  sum”  technique, 
for  sequential  input  gaussian  pressure  profiles,  with  peak  amplitudes  of  10  kPa,  and 
center  frequencies  of  1.8,  3.6  and  5.4  MHz,  respectively.  The  fourth  response  is  that 
obtained  using  Transducer  C  when  the  input  is  the  5.4  MHz  gaussian  pulse. 
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requires  only  the  use  of  simple  delays  similar  to  those  conventionally  used  in  an  imaging 
process.  The  “delay  and  sum”  function  may  be  implemented  in  either  hardware  or  software. 

4.4  Stability 

The  stability  of  the  “optimized  excitation”  method  for  transmission  is  considered.  Since 
the  technique  involves  modifying  more  than  one  input  in  a  manner  which  relies  on  the  ability  to 
model  the  practical  system,  it  is  natural  to  be  concerned  that  small  deviations  in  the  input 
quantities,  resulting  from  any  one  of  a  number  of  possible  causes,  will  result  in  a  relatively 
distorted  output. 

The  stability  of  the  technique  with  respect  to  deviations  in  the  applied  phase 
characteristics  has  been  examined.  Figure  31  illustrates  the  dependence  of  the  magnitude  of 
output  pressure  on  frequency  and  phase  angle,  9i  •  62.  for  Transducer  C.  The  complete  range 
of  phase  angles  has  been  considered.  A  smooth,  rounded  and  continuous  surface  is  evident. 
The  absence  of  discontinuities  indicates  that  small  phase  errors  will  not  result  in  substantial 
errors  in  the  output  quantity.  The  surface  can  be  made  more  even,  in  terms  of  magnitude,  by 
imfxoving  the  acoustic  match  between  the  piezoelectric  layers  and  the  backing  or  load  media. 

The  influence  of  errors  in  layer  dimension,  layer  material  properties,  layer  stimulation 
level  and  the  impact  of  a  finite  epoxy  bond  layer,  have  also  been  considered.  It  was  supposed 
that  Gaussian  pulses,  with  peak  amplitude  10  kPa,  and  center  frequencies  of  1.8  and  3.6  MHz, 
were  to  be  generated  using  Transducer  C.  These  frequencies  correspond  to  the  fundamental 
and  second  harmonic  of  the  combined  two  layer  system.  Hence,  the  3.6  MHz  pulse  required 
that  the  excitation  functions  to  layers  and  1  and  2  be  approximately  antiphase  with  respect  to 
each  other.  Figure  32,  illustrating  the  effect  of  reduced  layer  2  stimulation  applied  to  layer  2, 
indicates  mirumal  distortion  and  modest  attenuation  of  the  signal.  By  contrast.  Figure  33, 
illustrating  the  effect  of  reducing  the  thickness  of  layer  2  by  20%,  indicates  that  some  distortion 
has  been  introduced.  As  might  have  been  anticipated,  this  distortion  is  greater  at  the  higher 
frequency,  where  the  phase  characteristics  of  the  applied  signals  are  more  diverse.  Similarly, 
tlw  effect  of  altering  the  material  prc^rties  of  layer  2,  by  increasing  density  by  20%,  illustrated 
in  Figure  34,  indicates  more  distortion  at  the  higher  frequency.  Figure  35  illustrates  the 
consequence  of  the  most  probable  practical  defect.  The  introduction  of  a  5  p.m  epoxy  bond 
layer  between  layers  1  and  2  results  in  minor  distortion  at  the  higher  frequency. 
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Figure  31 

Surface  profile  indicating  dependence  of  pressure  output  magnitude  on  phase  angle, 
01  -  02,  and  frequency  for  Transducer  C, 
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Figure  32 

Theoietical  output  oi  Transducer  C  when  the  required  functions  are  sequential  1 .8 
and  3.6  MHz  pulses,  respectively.  The  desired  peak  amplitude  is  10  k^.  In  this 
case,  the  stimulation  of  layer  2  was  artificially  i^uced  by  20%.  The  response  is 
virtually  identical  to  the  desired  response  except  fn*  a  reduction  in  amplitude  of 
approximately  10%. 


Figure  33 

Theoietical  output  of  Transducer  C  when  the  requited  functions  are  sequential  1.8 
and  3.6  MHz  pulses,  tespecti\ely.  The  desired  peak  amplitude  is  10  kJ^.  In  this 
case,  the  thickness  of  layer  2  was  ariinciaUy  reduced  by  20%. 
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Rgure  34 

Theoretical  ouqiut  of  Transducer  C  when  the  required  functions  are  sequential  1 .8 
and  3.6  MHz  pulses,  reflectively.  The  desired  peak  amplitude  is  10  k)^.  In  this 
case,  the  density  of  layer  2  was  artificially  increased  by  20%. 


Figure  35 

Theoretical  output  ttf  Transducer  C  when  the  required  functions  are  sequential  1.8 
and  3.6  MHz  pulses,  respectively.  The  desired  peak  amplitude  is  10  kl^.  In  this 
case,  a  S  mm  epoxy  bond  layer  has  been  includ^  in  the  simulation. 
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The  deviatkxis  which  have  been  considered  here  are  relatively  extreme  and  the  resulting 
distortions  have  generally  been  minor.  It  can  be  concluded  that  the  technique  is  generally  stable 
but  that  the  imptntance  o(  quantifying  the  structure  accurately  is  more  signiticant  when  higher 
frequencies  are  employed. 

4.4  OonclusiQns  on  Multiple  Acttvc  Layg  Transducers 

The  sensitivity  and  bandwidth  characteristics  of  thickness  mode  transducers  can  be 
enhanced  by  using  multiple  active  layers.  The  variety  of  design  variables,  including 
dimensions,  material  properties  and  electrical  stimulation,  offer  the  designer  immense 
versatility.  It  is  impossible  to  determine  a  design  which  is  optimal  for  all  applications.  In 
transmission,  any  deficiencies  in  the  physical  design  may  be  compensated  for  by  adjusting  the 
stimulaticMi  functions.  However,  for  optimized  reception  characteristics  there  is  a  strong 
motivation  for  using  a  design  which  includes  a  backing  that  is  acoustically  matched  to  the 
piezoelectric  layers.  The  mathematical  technique  for  determining  the  excitation  functions  is 
straightfrvward  and  the  required  electronic  hardware  can  be  implemented  using  standard,  low 
cost  components.  The  technique  may  be  extended  to  pulse-echo  operation.  The  stability  of  the 
technique  has  been  examined  and  been  found  to  be  generally  good.  Although  only  two  layer 
cases  have  been  considered  here,  it  is  possible  to  extend  the  concept  by  using  more  active 
layers,  and  consequently  obtain  even  higher  bandwidths. 

5.  ENHANCED  1:3  COMPOSITE  TRANSDUCERS 

Composite  transducers,  comprising  aligned  piezoceramic  pillars  embedded  in  a 
polymer  matrix,  have  been  studied  extensively  and  have  found  applications  in  medical 
ultrasound,  non-destructive  examination  and  high  frequency  sonar.  However,  their 
efficimey  is  limited  by  the  positive  Poisson's  ratio  observed  in  all  conventional  polymers. 

It  has  recently  been  shown  that  device  perfcxmance  could  be  significantly  improved  if  a 
filler  possesang  a  negative  Poisson’s  ratio  could  be  used  [41].  Since  the  use  of  negative 
Poisson's  ratio  materials  is  not  readily  feasible  at  present,  a  compromise,  which  in  effect 
simulates  a  zero  Poisson’s  ratio  material,  is  presented  here.  Aligned  compressible 
inclusions  are  framed  in  a  standard  polymer  filler.  These  absorb  lateral  strain  without 
associated  Poisson's  ratio  cross-coupling  into  thickness  strain.  Typically,  these  inclusions 
are  gas  filled  holes  drilled  into  the  polymer.  The  resulting  structure  is  illustrated  in  Fig.  36. 
Notice  that  when  standard  notation  [42]  is  used,  the  resulting  composite  is  a  1:3:1 
composite  in  contrast  to  the  conventional  1:3  composite.  The  first  phase  (ceramic)  is 
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Vertically  aligned 
ceramic  pillars 


Figure  36 

New  composite  structure  with  aligned  holes. 


connected  only  in  the  thickness  direction.  The  second  phase  (polymer)  is  connected  in  three 
dimensitHis  and  third  phase  (gas)  is  connected  only  in  the  thickness  direction. 

5.1  Thickness  Mate  Tiansirnttcrs  and  Rwcivgrs 

In  a  conventitMial  1:3  composite  transducer,  comprising  ceramic  pillars  and  a 
polymer  matrix,  the  situation  illustrated  diagrammatically  in  Fig.  37  may  arise.  In 
transmission,  the  ceramic  is  die  active  displacen%nt  initiating  phase.  The  Poisson's  ratio  of 
the  ceramic  is  typically  +0.35  for  Lead  Zirconate  Titanate  (PZT)  and  the  polymer's 
Poisson's  ratio  varies  between  +0.35  for  epoxy  and  +0.47  for  polyurethane.  Hence, 
significant  antiphase  cross  coupled  strain  results  in  the  thickness  direction.  Consequently, 
the  net  thickness  displacement  of  the  polymer  may  be  antiphase  with  respect  to  that  of  the 
ceramic.  A  degree  of  cancellation  of  thickness  displacement  occurs  and  the  net  surface 
displacement  is  reduced.  This  has  had  a  significant  bearing  on  polymer  selection.  A 
compliant  polymer  possessing  a  low  Poisson’s  ratio  would  be  ideal.  Unfortunately, 
compliant  polymers  generally  have  high  Poisson’s  ratios.  For  this  reason,  less  compliant 
epoxy  is  generally  favored  in  preference  to  polyurethane,  which  is  more  compliant.  The 
properties  of  these  two  polymers  are  presented  in  the  Appendix  I.  Clearly,  the  extent  of  the 
problem  is  a  complex  function  of  geometry  and  material  properties.  Additionally  the  static 
consideration  presented  here  is  modified  in  the  high  frequency  thickness  resonant  case. 
Furthermore,  practical  transducers  are  three  dimensional  structures.  A  numerical  analysis, 
namely  finite  element  analysis,  is  required  to  obtain  a  useful  model  of  the  phenomena 
present. 

Nevertheless,  if  we  consider  the  modified  transducer,  illustrated  diagrammatically 
in  Fig.  38,  with  compressible  gas  filled  inclusions,  it  is  immediately  evident  that  all  induced 
lateral  strains  in  the  polymer  are  absorbed  by  the  inclusions.  Immediately,  a  result  is 
obtained  which  simultaneously  leads  to  improved  thickness  displacement  uniformity  and 
has  the  potential  for  increased  electromechanical  coupling  efficiency.  The  only  coupled 
thickness  motion  in  the  polymer  is  due  to  thickness  extensional  coupling  across  the 
ceramic/  polymer  boundary.  Naturally,  this  thickness  extensional  coupling  in  the  polymer 
is  in  phase  widi  the  ceramic.  Considerable  importance  has  been  attached  to  ensuring  that  the 
pacing  between  adjacent  pillars  is  not  sufficiently  great  that  parasitic  shear  standing  waves 
can  resonate  in  the  vicinity  of  the  thickness  operating  frequency  [43,44].  This  is  a  different 
phenomenon,  occuring  only  in  the  dynamic  regime.  Notice  that  the  presence  of  the  holes 
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Exaggerated  illustration  of  anti-phase  motion  resulting  from  positive  Poisson’s  ratio 
in  ceramic  and  polymer  phases. 
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Figiro  38 

Exagiprated  illustration  of  cmnpressible  inclusions  absorbing  strains  without  cross- 
coupling  to  thkki^  direction. 
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will  probably  have  a  disruptive  effect  on  the  lateral  standing  waves  normally  occurring  in  a 
conventional  composite  transducer. 

A  very  important  associated  effect  is  that,  since  the  polymer  resists  lateral 
displacement  of  the  ceramic  (Fig.  37),  a  significant  degree  of  lateral  clamping  pressure  is 
applied  to  the  lateral  surfaces  of  the  ceramic.  This  has  the  effect  of  reducing  the  potential 
electrcmiechanical  coefficient  (k33)  of  a  PZT  5A  [45]  bar  from  0.70  towards  the  “laterally 
clamped”  coefficient  (kt),  which  is  0.49.  This  decrease  will  be  reduced  in  the  new  1:3:1 
structure.  A  less  significant  associated  effect  is  that  as  lateral  stiffness  is  reduced  in  the 
composite  so  also  is  stiffness  in  the  thickness  direction.  Hence  a  modest  reduction  in  the 
thickness  nKxle  rescmant  frequency  is  obtained.  This  implies  that  the  acoustic  impedance  of 
the  composite  is  also  reduced  in  the  new  structure,  lliis  is  a  desirable  effect  for  the  majority 
of  applications. 

In  the  current  program  of  work,  finite  element  analysis  has  been  performed  using 
the  commercially  available  ANSYS  code  [46].  A  full  discussion  of  the  development  of  the 
thecffy  relating  to  the  finite  element  analysis  of  piezoelectric  materials  may  be  found  in  the 
literature  [47,48].  Furthermore,  the  application  of  the  ANSYS  finite  element  code  to 
ccmqxrrite  transducers,  with  experimental  validation,  is  discussed  in  detail  in  Reference  49. 
In  all  cases,  unless  specifically  noted  otherwise,  the  natural  resonant  frequencies,  and  the 
nature  of  the  associated  modal  displacements,  have  been  evaluated  using  eigenvalue 
analysis  [47].  Relative  displacement,  as  opposed  to  absolute  displacement,  information  is 
genoated.  Repeating  the  aiudysis  for  both  short  circuit  and  open  circuit  electrode  conditions 
permits  the  resonant  and  antiresonant  frequencies  to  be  evaluated  [49].  Subsequently,  the 
electromechanical  coupling  coefficient  associated  with  a  particular  tiKxle  may  be  evaluated. 
An  additional  indication  of  the  practical  significance  of  a  particular  mode  is  obtained  by 
considering  the  relative  thickness  displacements  of  the  nodes  on  the  transducer  surface. 
Qearly  a  mode  in  which  all  nodes  displace  by  a  similar  extent  is  of  more  practical 
usefulness  than  a  traxie  in  which  certain  tiKxies  displace  in  one  direction  and  others  in  the 
qpposite  direction. 

Since  the  composite  structures  considered  here  have  quaner  symmetry  in  the  X-Y 
(latoal  plane)  and  strongly  coupled  trxxfes  are  symmetric  in  the  Z  (thickness  direction),  it  is 
only  r»cessary  to  model  a  (me  eighth  section  of  the  unit  cell.  Naturally,  care  is  taken  to 
inqwse  the  approbate  boundary  conditions.  No  perpendicular  displacement  of  nodes  lying 
(m  a  symmetry  plane,  with  respect  K>  the  plane,  is  permined.  The  material  properties  of  the 
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PZT  5A  ceramic,  a  typical  epoxy  polymer,  and  polyurethane  are  presented  in  Appendix  I. 
The  performance  improvement  induced  by  the  modification  discussed  in  this  paper  is  far 
rntne  apparent  when  polyurethane  is  employed.  Consequently,  the  new  design  is  most 
tqiplicable  to  devices  incorporating  this  material.  Comparisons  with  results  obtained  when 
epoxy  is  used  are  made,  as  required. 

Figure  39  illustrates  a  one  eighth  section  modal  displacement  representation  of  a 
25%  ceramic  volume  fraction  composite  transducer.  Plot  A,  referring  to  a  device  containing 
a  conventional  polymer  possessing  a  Poisson's  ratio  of  +0.47,  indicates  uneven  thickness 
displacement,  with  a  degree  of  antiphase  motion.  By  contrast.  Plot  B  indicates  the  results 
for  an  identical  device  except  that  the  Poisson's  ratio  of  the  polymer  is  0.0.  In  accordance 
with  the  argument  set  out  above,  it  is  observed  that  in  case  B  a  more  even  displacement 
occurs.  The  immediate  aim  is  to  experimentally  realise  an  approximation  to  this  case  by 
introducing  compressible  inclusions  into  the  polymer.  A  practically  viable  method  would  be 
to  drill,  or  otherwise  form,  holes  in  the  material.  Alternatively,  one  could  saw  along  the  X- 
Y  directions  between  adjacent  ceramic/  polymer  cells.  This  would  decouple  each  ceramic/ 
polymer  cell  from  neighboring  cells  and  virtually  zero  net  lateral  stress  conditions  would 
exist.  Three  hole  geometries  are  illustrated  in  Fig.  40  in  addition  to  the  inter  pillar  sliced 
geometry.  For  convenience  these  geometries  will  subsequently  be  refetred  to  as  “diagonal,” 
“parallel,”  “combined,”  and  “sliced” — as  defined  in  the  figure.  It  is  immediately  evident 
that  a  virtually  infinite  range  of  ceramic/polym^/hole  geometries  is  possible.  A  limited 
range  is  considered  here  with  a  view  to  obtaining  a  limited  understanding  of  the  phenomena 
present  so  that  informed  design  decisions  can  be  made. 

A  substantial  simplification  of  the  finite  element  mesh  is  afforded  if  it  is  assumed 
that  both  the  ceramic  pillars  and  the  holes  are  rectangular  in  cross-section.  This  avoids  the 
more  complex  geometry  associated  with  frtting  circular  features  within  rectangular  frames. 
In  practice  the  ceramic  pillars  are  normally  square  but  the  holes  would  probably  be  circular. 
The  influence  of  the  exact  shape  is  of  little  significance  in  the  examples  presented  here  and 
the  added  complexity  of  a  more  realistic  mesh  is  not  justifred.  Analysis  of  the  relative 
performance  of  circular  and  square  section  ceramic  pillar  composites  described  in  Reference 
49  serves  to  justify  this  approximation. 

One  might  suppose  that  an  alternative  approach  would  be  to  use  a  ceramic  with  a 
low  Poisson's  ratio.  Unfortunately,  modified  Lead  Titanate  (PT),  the  only  ceramic  with 
such  a  characteristic,  has  a  relatively  low  electromechanical  coupling  coefficient  in  the  bar 
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Modal  displacement  at  the  thickness  resonant  frequency  in  a  25%  ceramic  volume 
fraction  composite  containing  two  types  of  polymer.  A  one-eighth  section  is 
illustrated,  llie  ceramic  pillar  is  0.5  mm  wide  and  2  mm  high. 
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(a)  Polyurethane  polymer  (Poisson’s  ratio  =  0.47) 

(b)  Mo^fied  polymer  (Poisson’s  ratio  =  0.0) 
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Cross-sections  of  various  hole  geometries. 
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mode  In  PT,  k33  =  0.50,  whereas,  for  PZT  k33  =  0.70.  Additionally,  PT  has  a  low 
dielectric  permittivity.  This  problem  is  compounded  in  a  composite  since  the  net 
permittivity  is  further  reduced  due  to  the  replacement  of  part  of  the  ceramic  with  polymer  in 
the  structure. 

Results  obtained  using  tinite  element  analysis  were  obtained  for  ceramic  volume 
fractions  of  both  25%  and  50%.  In  order  to  shift  the  effect  of  parasitic  lateral  modes  to 
frequencies  well  separated  from  the  thickness  resonant  frequency,  a  pillar  width  to  height 
(aspect)  ratio  of  0.1  has  been  used. 

In  the  first  theoretical  test  case,  the  effect  produced  by  increasing  the  number  of 
holes  in  the  structure  is  investigated.  The  same  basic  mesh  geometry  is  used  so  that  relative 
comparisons  of  nodal  displacements  are  valid.  In  practice  the  locations  of  the  nodes  were 
constant  but  an  increasing  number  of  polymer  elements  were  replaced  by  empty  spaces 
simulating  holes.  In  the  25%  ceramic  volume  fraction  case,  each  hole  occupied  6.25%  of 
the  total  volume,  whereas,  in  the  50%  ceramic  volume  fraction  case,  each  hole  occupied 
2.14%.  These  geometries  were  chosen  so  as  to  provide  a  regularly  spaced  geometric  mesh. 

Figure  41  illustrates  the  influence  of  hole  geometry  on  the  thickness  mode 
electromechanical  coupling  coefficient.  For  both  the  25%  and  50%  ceramic  volume  fiction 
composites  a  significant  improvement  is  achieved  simply  by  inserting  one  hole  diagonally 
between  adjacent  pillars.  Smaller  improvements  are  obtained  with  increasing  numbers  of 
holes.  As  would  be  expected,  the  “sliced”  geometry,  which  is  totally  unconstrained  in  the 
lateral  direction,  provides  the  optimal  result. 

A  measure  of  the  degree  of  nodal  displacement  uniformity  is  presented  in  Fig.  42. 
The  standard  deviation  of  perpendicular  nodal  displacement  on  the  transducer's  vibrating 
surface  has  been  calculated.  Since  modal  analysis  generates  only  relative,  not  absolute, 
displacement  information,  the  largest  displacement  was  normalized  to  1.0  and  all  other 
displacements  were  scaled  by  a  similar  factor.  Subsequently,  the  statistical  standard 
deviation  of  these  normalized  displacements  was  calculated.  Although  there  is  no 
monottmic  improvement  in  the  uniformity  of  nodal  displacement,  it  is  readily  {^parent  that 
a  significant  improvement  is  obtained  by  fabricating  one  or  more  holes.  The  exact  nature  of 
the  nodal  displacement  is  a  complex  function  of  geometry  and  consequently  it  is  not 
possible  to  predict  intuitively  that  the  uniformity  of  nodal  displacement  will  always  improve 
with  increasing  numbers  of  holes.  It  is  evident  that  the  “sliced”  structure,  which  is 
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Electromechanical  coupling  coeff. 


Figure  41 

Theoretical  dependence  of  thickness  mode  coupling  coefficient  on  hole  geometry 

-  25%  Ceramic  volume  faction 

.  50%  Ceramic  volume  fraction 
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Standard  deviation  of  nodal  disp 


Figure  42 

Theoretical  dependence  of  nodal  displacement  uniformity  (thickness  direction)  on 
hole  getxnetry 

— -  25%  Cteramic  volume  fraction 

.  50%  Ceramic  volunnc  fraction 
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c(»iq>ietely  unconstrained  on  the  lateral  surfaces  of  each  pillar/polymer  unit  cell,  does  not 
always  vibrate  uniformly.  However,  the  high  electromechanical  coupling  coefficient 
associated  with  this  structure  compensates  to  a  lar^  extent  for  this  potential  deflcioicy. 

Figure  43,  illustrating  Icmgitudinal  wave  velocity  in  the  thickness  direction,  as  a 
fiinctiiMi  of  hole  gecxnetiy,  indicates  a  general  reduction  in  vdocity  as  the  number  of  holes 
is  increased.  This  parameter  was  calculated  from  the  anti-rescmant  frequency  using  the 
assumption  that  this  corresponds  to  a  V2  resonant  condition.  This  characteristic  is  in 
accordance  with  intuitive  expectation.  Increasing  numbers  of  holes  result  in  reduced  lateral 
'clanqnng'  which  cross  couples  to  the  thickness  direction  in  the  form  of  reduced  stiffness 
in  the  thickness  direction,  as  noted  previously. 

Up  to  this  point  the  analyses  have  considered  one  mesh  for  the  25%  ceramic 
volume  fraction  device  and  another  for  the  50%  ceramic  volume  fraction  device.  Hence,  as 
the  number  of  holes  increases,  so  also  does  the  void  volume  fraction.  It  is  instructive  to 
consider  independently  the  effects  of  increasing  numbers  of  holes,  and  increasing  void 
fraction. 

Figure  44  illustrates  the  dependence  of  electromechanical  coupling  coefficient  on 
void  volume  fraction  for  the  25%  ceramic  volume  fraction  composite.  Analysis  was 
performed  for  the  diagonal,  parallel  and  combined  hole  geometries.  The  coupling 
coefficient  associated  with  a  conventional  structure  and  with  a  sliced  structure  are  also 
illustrated.  It  is  evident  that  the  coupling  coeffrcient  is  enhanced  as  the  void  fraction 
increases.  However,  the  improvement  levels  off  at  higher  void  fractions.  It  is  also  evident 
that  the  diagonal  hole  design  is  the  most  effective  of  the  three  geometries.  This  may  be  due 
to  the  fact  that  the  hole  dimension  is  greatest  when  only  a  single  hole  is  employed.  The 
results  fOT  the  50%  ceramic  volume  fraction  composite,  illustrated  in  Fig.  45  demonstrates 
similar  trends.  It  is  physically  impossible  to  achieve  a  large  void  fraction  with  a  single 
square  hde  in  a  50%  ceramic  volume  fraction  composite.  Hence,  some  corresponding  data 
points  are  missing.  It  is  apparent  that  substantial  improvements  are  achieved  with  smaller 
void  fractions  in  the  case  of  the  50%  ceramic  volume  fraction  composites.  This  is  probably 
related  to  the  fact  that  the  polymer  volume  fraction  is  also  less  in  this  case. 

Similar  trends  were  observed  when  the  analysis  was  performed  for  epoxy  filled 
conqxisites.  The  only  difference  was  that  the  degree  of  enhancement  of  thickness  mode 
coupling  coeffrcient  was  approximately  one  third  of  that  shown  in  Figs.  41, 44  and  45. 
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Figure  43 

Theoretical  dependence  of  longitudinal  wave  velocity  (thickness  direction)  on  hole 
geometiy 

-  25%  Ceramic  vdume  fraction 

.  50%  Ceramic  volume  fraction 
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Electromechanical  coupling  coeff. 


Figure  44 

Theoretical  dependence  of  thickness  electromechanical  coupling  coefficient  on  void 
volume  fraction  in  a  25%  ceramic  volume  fracticm  composite 
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Figure  45 

Theoretical  dependence  of  thickness  electromechanical  coupling  coefficient  on  void 
volume  fraction  in  a  50%  ceramic  volume  fraction  composite 
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Two  prototype  enhanced  composite  transducers  have  been  manufactured  and  tested. 
The  geometries  of  these  are  described  in  Appendix  ni.  Unfortunately,  only  epoxy  was 
available  as  a  polymer  and  hence  any  improvements  demonstrated  are  less  than  might  have 
been  obtained  had  polyurethane  been  used.  Also,  practical  limitations  on  the  size  of  the 
holes  constrained  the  overall  geometry.  Since  the  holes  were  hand-drilled,  it  was  necessary 
to  use  relatively  wide  spacings  between  the  pillars.  This  gave  rise  to  the  presence  of 
atkiitional  resonant  modes  slightly  higher  than  the  thickness  resonant  frequency.  This 
prevented  accurate  measurement  of  the  electromechanical  coupling  coefHcients.  The 
transmissicm  performance  was  measured  with  a  hydrophone  placed  on  axis,  in  the  far  Held, 
in  a  water  tank.  A  sinusoidal  burst,  at  the  center  frequency  and  approximately  ten  cycles 
long,  was  applied  to  each  transducer.  The  responses  were  normalized  with  respect  to  the 
response  of  the  conventicxial  composite  transducer  and  are  indicated  in  Table  I. 

TABLE  I 

Measured,  normalized  transmission  and  reception 
performance  of  Transducers  I  and  II 


Hole  Configuration 

Conventional 

Diagonal 

Parallel 

C^ombined 

Transmissicxi  (Txdr  A) 

1.00 

1.11 

1.11 

1.15 

Transmission  (Txdr  B) 

1.00 

1.18 

1.18 

1.18 

Reception  (Txdr  A) 

1.00 

1.25 

1.31 

1.25 

Reception  (Txdr  B) 

1.00 

1.12 

1.12 

1.20 

Improvements  in  sensitivity,  in  the  range  10%  to  20%  were  achieved  in  these  cases.  The 
reception  response  was  measured  using  a  wideband  transmitting  transducer  operating  in  the 
far  field,  with  the  distance  between  the  transmitter  and  the  detecting  transducer  greater  than 
the  sum  of  the  near  field/  far  field  transition  ranges  of  both  transducers.  A  fixed  amplitude, 
sinusoidal  burst,  at  the  center  fi^uency,  was  applied  to  the  transmitter.  The  response  of 
each  receiving  transducer  was  monitored  and  the  amplitude  normalized  with  respect  to  the 
response  of  the  conventional  design.  Improvements  in  the  range  10%  to  30%  were 
obtained  by  using  the  new  designs.  The  results  clearly  indicate  that  improved  performance 
is  obtainable.  However,  it  will  be  desirable  to  develop  manufacturing  techniques  which 
would  permit  accurate  manufacture  of  pillar/polymer/hole  geometries  finer  than  those 
available  at  this  time.  It  is  also  evident  that,  for  the  purposes  of  demonstrating  the 
effectiveness  of  the  technique,  it  would  be  desirable  to  use  polyurethane  rather  than  epoxy. 
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In  conclusion,  improvement  in  composite  transducer  performance,  in  terms  of 
elecmnnechanical  coupling  coefficient,  surface  displacement  uniformity  and  reduced 
longitudinal  wave  velocity,  is  achievable  by  employing  aligned  compressible  inclusions.  It 
would  appear  that  this  improvement  is  significant  when  polyurethane  is  employed. 
However,  the  improvement  observed  when  epoxy  filled  composites  were  considered  was 
far  less.  This  is  explained  by  the  observation  that  Polyurethane  possesses  a  far  higher 
Poisson's  ratio.  Nevertheless,  there  has  been  a  long-standing  interest  in  soft  polymers  for 
use  in  flexible  conq)osite  structures.  These  soft  polymo's  have  been  largely  abandoned  as  a 
component  in  ccnnposites  because  of  their  large  Poisson's  ratio.  The  new  technique 
provides  a  means  for  circumventing  the  problem  of  a  high  Poisson's  ratio. 

5.2.  Hydrophones 

The  physical  design  of  a  hydrophone  using  the  new  composite  is  identical  to  that 
discussed  above  for  the  thickness  mode  transmitter/  receiver.  However,  the  principle  of 
operation  is  quite  different.  In  the  hydrostatic  mode,  the  net  charge  genoated  by  hydrostatic 
pressure  is  (tetermined  by  dh  charge  coefficient  This  parameter  is  related  to  d33  and  d3 1  in 
the  following  maimer. 

dh*d33  +  2d3i  (27) 

The  hydrostatic  charge  coefficient  is  dh  (units  C7N  or  m/V),  while  d33  is  the 
coefficient  relating  thickness  induced  charge  to  applied  thickness  stress  and  d3i  is  the 
coefficient  relating  thickness  induced  charge  to  applied  latoal  stress.  PZT  materials  are  very 
efficient  charge  generators.  However,  for  PZT  5A,  d33  =  374  x  10*12  m/V  and  d3i  = 
-171  X  10*12  m/V.  Therefore,  the  dh  coefficient  is  only  32  x  10*12  m/V.  It  is  evident  that 
the  objective  in  the  design  of  hydrophones  incorporating  PZT  should  be  to  decouple  the 
thickness  and  lateral  stresses.  A  composite  design  achieves  this  to  an  extent  since  the 
polymer  phase  may  serve  to  absorb  lateral  strains.  However,  if  the  polymer  has  a  high 
Poisson's  ratio,  thickness  strain  is  still  strongly  coupled  to  lateral  strain.  The  phase 
relatitmship  of  this  induced  strain,  with  respect  to  the  thickness  strain,  implies  that  the  net 
induced  charge  is  reduced.  One  technique  for  reducing  the  extent  of  cross  coupling  of  strain 
in  the  polymer  phase  is  to  introduce  porosity  [50].  A  foam  type  structure  achieves  this 
effect  successfully  but  it  is  vulnerable  to  collapse  when  subjected  to  high  hydrostatic 
pressure  [50].  Another  approach  involves  weaving  glass  fibers  into  the  composite  matrix. 
These  absorb  smne  of  the  lateral  stress  but  enhance  the  net  charge  generating  effect  to  a 
lesser  extent  dian  the  foamed  polymer  techniqiK.  Another  design  uses  a  substantial  air  void 
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sunounding  all  the  ceramic  pillars  [SI].  High  sensitivity  is  achieved  but  it  is  evident  that  the 
design  is  vulnerable  to  damage  at  high  hydrostatic  pressures. 

The  new  hydrt^hone  design  comprises  of  the  same  aligned  ccHnpressible  inclusions 
as  discussed  in  the  previous  section.  Perhaps  the  greatest  attribute  of  the  technique  is  that  a 
highly  anisotix^ic  net  polymer/  void  characteristic  is  achieved.  Pressure  waves  impinging 
perpendicular  to  the  thickness  facing  surface  encounter  a  “parallel  springs"  structure  of 
alternating  ceramic,  polymer  and  gas.  The  ceramic  is  the  dominant  “spring”  and  hence  it 
preferentially  absmbs  the  induced  stress.  However,  in  the  lateral  direction,  an  incoming 
wave  encounters  a  “series  spring”  structure  composed  of  ceramic,  polymer  and  gas. 
Naturally,  the  gas  dominates  in  the  series  case  and  absorbs  virtually  all  induced  str  n. 
Hence,  optimal  use  is  made  of  a  highly  anisotropic  composite  structure. 

The  effect  of  introducing  holes  in  a  composite  has  been  evaluated  theoretically  using 
finite  elonent  analysis.  One  basic  design  has  been  tested.  The  ceramic  volume  fraction  was 
25%  and  the  ceramic  pillars  were  4  mm  long  and  1  mm  across.  Devices  containing  epoxy 
and  polyurethane  have  been  evaluated  separately.  Since,  in  a  practical  system,  it  is 
necessary  to  have  a  layer  over  the  surface  to  protect  the  electrode  and  to  prevent  water 
ingress  into  the  holes,  the  design  has  also  been  tested  with  a  0.5  mm  thick  layer  on  the  top 
and  bottom  surfaces.  The  holes,  where  they  have  been  introduced,  were  0.5  mm  across 
and  extoided  throughout  the  thickness  of  the  active  composite  material.  It  is  anticipated  that 
the  design  will  be  resistant  to  high  hydrostatic  pressure  since  the  holes  need  only  possess  a 
small  cross-sectional  area.  Additionally,  unlike  the  foamed  structures  considered 
previously,  these  are  txdered  structures  and  hence  weak  points  can  be  avoided  by  judicious 
design.  The  finite  element  analysis  was  performed  by  applying  a  known  static  pressure  in 
three  planes.  The  charge  induced  was  related  to  the  applied  force  to  yield  the  effective  dh 
coefficient  The  technique  was  initially  tested  using  a  reference  ceramic  sample.  The  d33 
coefficient  calculated  using  this  technique  matched  the  manufucturer's  tabulated  value  [45]. 

The  theoretical  results  obtained  are  presented  in  Table  U.  It  may  be  observed  that 
very  substantial  increases  in  the  dh  coefficient  are  obtainable  Generally,  the  two  hole  design 
A)es  not  produce  results  substantially  different  from  the  one  hole  design.  However,  using 
the  combirted  hole  design  yields  better  results.  The  “sliced”  design  produces  by  far  the  best 
results.  In  this  case,  the  lateral  force  is  entirely  decoupled.  This  design  would  be 
sutetantially  more  vulnoable  to  collapse  at  high  hydrostatic  pressure.  It  is  also  evident  that, 
as  in  the  case  of  die  thickness  resonant  transducers,  a  far  greater  enhancement  is  obtained 
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Depenctence  of  ri 


TABLE  II 

ical  dh  coefficient  on  composite  geometry  and  material, 
(dh  values  in  units  IQ-^^  C/N) 


No  layer  With  layer 


Polyurethane 

Epoxy 

Polyurethane 

Epoxy 

No  hole 

1.87 

8.05 

4.09 

9.07 

Diagcxial 

13.27 

19.12 

17.40 

20.76 

Parallel 

10.51 

16.67 

19.45 

22.07 

Ccxnbined 

27.69 

32.24 

38.73 

38.59 

Sliced 

272.97 

248.06 

288.17 

259.34 

by  using  holes  when  the  polymer  is  polyurethane  rather  than  epoxy.  This  is  because 
polyurethane  is  largely  incompressible  and  hence  strains  are  strongly  coupled  between  the 
lateral  and  thickness  directions  when  the  polymer  is  in  bulk  form.  At  least  in  the  case  of 
epoxy,  the  presence  of  a  protective  layer  does  not  significantly  modify  the  degree  of 
improvement  that  tme  might  expect  from  using  the  technique.  When  it  is  considered  that  the 
commonly  used  dhgh  coefficient  [52]  is  related  to  dh^,  it  is  clearly  evident  that  the 
technique  presented  here  is  capable  of  producing  very  substantial  improvements  in 
performance. 

6.  PROGRAM  OF  COLLABORATION  WITH  FMI 

A  collaboration  program  with  FMI,  Maine,  was  established.  This  provided 
StanfOTd  with  a  facility  for  obtaining  samples  and  experimental  results,  and  FMI  with  a 
means  of  improving  their  understanding  of  composite  behavior  with  a  view  to  establishing 
optimal  practical  designs.  FMI  has  a  manufacturing  technology  with  greater  flexibility  than 
that  cf  the  conventional  dice-and-  fill  technique.  Specifically,  elements  can  be  inserted  in  a 
polymer  substrate  in  a  completely  arbitrary  two-dimensional  pattern. 

Various  phenomena  were  to  be  investigated; 

1 .  Reducing  parasitic  mode  behavior  using  the  following  techniques. 

a)  Specialized  geometries  resulting  in  canceling  of  piezoelectric  coupling  to 
specific  bandedge  resonances. 

b)  Statistically  varying  inter-pillar  dimensions 

c)  Filler  ma^als 
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2.  Perfonnance  as  a  function  of  polymer  properties 

3.  Performance  of  new  “enhanced”  composite  with  aligned,  compressible,  gas-filled 

inclusions. 

The  starting  point  of  the  program  was  the  establishment  of  a  suitable  baseline 
design  possessing  well-characterized  parasitic  resonant  modes.  Sample  #16822,  using  a 
square  lattice,  exhibited  three  clearly  identiHable  resonances.  This  sample  has  a  ceramic 
volume  fraction  of  10%,  and  the  PZT-5H  pillars  are  1  mm  across.  FMI  HD-68  epoxy  was 
employed  as  the  filler.  The  sample  measures  38  x  38  mm  and  is  6.35  mm  thick.  The  three 
resonances  illustrated  in  Fig.  46,  at  progressively  higher  frequencies,  may  be  related  to  the 
parallel  lateral  mode,  the  thickness  mode,  and  the  diagonal  lateral  mode,  respectively. 
When  the  structure  was  analyzed  using  frnite  element  analysis,  the  corresponding  modes, 
illustrated  in  Fig.  47  were  obtained.  The  correlation  between  Figs.  46  and  47  was 
established  by  comparing  the  frequencies.  Note  that  the  illustrations  in  Fig.  47  are  for  a 
one-eighth  section  of  a  unit  cell,  with  the  centerline  of  the  cell  at  the  upper  comer  and  the 
horizontal  mid-plane  at  the  bottom. 

6.1  Reducing  parasitic  mode  activity 
a)  Specialized  Geometries 

The  first  modified  structure  involves  using  a  rectangular  rather  than  a  square  lattice. 
This  lattice  is  illustrated  in  Fig.  48.  The  principle  is  that  at  the  lowest  frequency  bandedge 
resonance,  there  is  V2  spacing  along  the  horizontal  axis  and  A/4  spacing  along  the  vertical 
axis.  Alternate  vertical  columns  vibrate  anti-phase  and  hence  the  integrated  piezoelectric 
coupling  is  zero.  At  the  next  higher  resonance,  with  A,  spacing  horizontally  and  A/2  spacing 
vertically,  alternate  horizontal  rows  vibrate  anti-phase.  Figure  49  illustrates  the  impedance 
plot  of  sample  #16714,  which  possesses  this  geometry.  It  is  clear  that  this  technique  does 
in  f«;t  result  in  reduced  lateral  mode  activity.  However,  some  remnant  resonant  activity  is 
still  present 

A  secoml  nxxlifred  structure  is  similar  to  that  illustrated  in  Fig.  48  except  that  pillars 
are  placed  in  the  centers  of  the  recu  gular  lattices.  Scattering  from  the  edges  of  the 
rectangular  lattices  are  1 80<’  out  of  phase  with  scattering  from  the  center  and  hence  the 
lowest  Bragg  order  is  suppressed.  Figure  SO  illustrates  the  plot  characteristic  of  sample 
#16715,  which  possesses  this  structure.  The  diagonal  resonance  has  been  completely 
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PN  16822 

A:  1Zl  B:  e  o  MKR  255  250.000  Hz 

A  MAX  10.00  Kn  MAG  1.39133  KQ 

B  MAX  80.00  deg  PHASE  38.1730  deg 


A  MIN  200.0  Q  START  100  000.000  Hz 

B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Figure  46 

Impedance  characteristic  of  the  baseline  sample  #16822  (Square  lattice) 
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(a) 


(b) 


Figure  47 

Finite  element  derived  modes  in  the  baseline  sample 
Modal  displacement  shape - 

Static  shape - 

a)  Parallel  mode 

b)  Thickness  mode 

c)  Diagonal  mode 
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PZT  COMPOSITE  ARRAY  DESIGN  #1 
PART  SIZE:  1.5'  X  1.5' 


PART  NUMBER  16714 


Figure  48 

Rectangular  lattice  design 


84 


PN  16714 

A:  iZi  B:  8  o  MKR  228  250.000  Hz 

A  MAX  10.00  KQ  MAG  1.25911  KQ 

B  MAX  80.00  deg  PHASE  46.6417  deg 


A  MIN  200.0  Q  START  100  000.000  Hz 

B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Figure  49 

Impedance  characteristic  of  sample  #16714  (Rectangular  lattice) 


PN  16715 

A:  IZl  B:  9  o  MKR  298  000.000  Hz 

A  MAX  10.00  Kfi  MAG  1.70322  KQ 

B  MAX  80.00  deg  PHASE  32.2099  deg 


A  MIN  200.0  Q  START  100  000.000  Hz 

B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Figure  50 

Impedance  characmistic  of  sample  #16715  (ModiOed  lattice) 
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eliminated.  The  lower  frequency,  parallel,  mode  remains.  The  phenomena  behind  this 
result  are  not  fully  understood  at  present. 

b)  Statistically  Varying  Geometries 

It  has  been  established  that  by  using  a  statistical  distribution  of  inter-pillar 
dimensions,  it  is  possible  to  eliminate  strong  parasitic  lateral  resonant  activity  occurring  at 
discrete  frequencies  [53].  The  previous  work  conclusively  proved  that  the  use  of  statistical 
distributions,  cerated  by  dicing  and  Hlling  with  quasirandomly  spread  dicing  cuts,  is 
effective.  The  two-dimensional  modeling  method  discussed  in  Section  7  suggests  that 
quasirandcxn  patterning  of  the  ceramic  element  locations  in  two  dimensions  should  be  more 
effective  in  mode  suppression.  This  kind  of  pattern,  not  realizable  by  the  dice  and  fill 
technique,  can  be  manufactured  using  FMI  technology.  Such  composites,  derived  ffom  the 
baseline  (#16822),  were  manufactured  by  FMI  with  variances  of  10%,  30%  and  50% 
applied  to  the  inter  pillar  spacings.  The  impedance  spectra  of  these  devices  are  illustrated  in 
Figs.  51, 52  and  53  respectively.  These  results  indicate  that  variances  of  10%  and  30%  are 
ineffective  in  terms  of  reducing  parasitic  resonant  activity.  However,  the  sample  with  50% 
variance  does  indicate  a  reduction  in  resonant  activity.  Previous  work  has  conclusively 
proven  that  the  use  of  statistical  distributions  is  effective  [53].  It  is  suspected  that  the  10% 
and  30%  results  did  not  appear  as  effective  as  they  might  otherwise  have  been  because 
tolerances  maintained  in  the  manufacturing  technique  currently  used  by  FMI  introduce  a 
degree  of  variance  into  the  basis  structure  as  a  matter  of  course.  The  degree  of  variance  has 
not  been  precisely  quantified,  since  it  is  operator  dependent.  In  future,  the  degree  of 
variance  may  be  more  tightly  controlled  by  the  use  of  complete  automation.  Nevertheless,  it 
would  appear  that  a  degree  of  variance  is  desirable  in  any  case. 

c)  Filler  Materials 

A  sample  similar  to  the  baseline,  (#16822)  except  that  it  included  glass  fiber 
reinfcHcement,  was  manufactured  and  tested.  Hber  reinforcement  has  a  beneficial  effect  chi 
the  hydrostatic  performance  of  composite  transducers  since  the  glass  absorbs  lateral 
stresses  thus  relieving  the  active  ceramic  pillars.  In  the  hydrostatic  mode  it  is  highly 
desirable  to  decouple  thickness  and  lateral  stresses  since  the  induced  charges  due  to  these 
separate  stresses  tend  to  cancel  each  other  out.  Additicmally,  reinforcement  may  lead  to 
inqnroved  mechanical  stability  and  result  in  attenuaticm  of  the  lateral  mode.  The  impedance 
spectrum  of  the  reinforced  device  is  illustrated  in  Fig.  54.  The  resonant  characteristics 
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16981. 

A:  IZl  B:  e  o  MKR  257  500.000  Hz 

A  MAX  ,10.00  KQ  MAG  1.6B063  KQ 

B  MAX  80.00  deg  PHASE _ 47 . 5274  deg 


A  MIN  200.0  n  START  100  000.000  Hz 
B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Hgure  51 

Impedance  characteristic  of  sample  #1^81  (10%  variance  applied  to  inter-pillar 
dimrosion) 
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B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Hgure  52 

cluBflcieristic  of  ssn^lc  #1^83  (30%  variance  applied  to  inter-pillar 

diniension) 
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FN  16986 

A:  IZl  B:  5  o  MKR  237  250.000  Hz 

A  MAX  10.00  Kfl  MAG  733.470  Q 

B  MAX  60.00  deg  PHASE _ -22.B519  deg 


M 

a 


A  MIN  500.0  Q  START  100  000.000  Hz 
B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Figure  53 

Impedance  characteristic  sample  #16986  (50%  variance  applied  to  inier-pillar 
dimeiision) 
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PN  16987 

A:  IZI  Q:  e  O  MKR  241  750.000  Hz 

A  MAX  10.00  Kfl  MAG  504.930  Q 

0  MAX  60.00  deg  PHASE  -17.7496  deg 


A  MIN  200.0  0  START  100  000.000  Hz 
B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Figure  54 

Impedance  characteristic  of  sample  #1^87  (Glass  fiber  reinforced  san^le) 
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differ  slightly  from  those  of  the  baseline  sample  (Fig.  46).  This  may  be  due  in  part  to  the 
fact  that  the  volume  fraction  of  leinfOTcement  fibers  was  only  5%  (#16987). 


A  second  composite  possessing  a  reinforcement  volume  fraction  of  10%  (#17491) 
was  also  manufactured.  In  this  case  Ceramaform  reinforcement  was  employed.  The 
impedance  characteristic  for  this  device,  illustrated  in  Fig.  55,  also  indicates  strong  parasitic 
resonant  behavior.  It  is  probable  that  the  resonant  behavior  would  be  more  significantly 
modified  if  higher  reinforcement  loading  fractions  were  used.  Given  the  vast  range  of 
possible  reinforcement  geometries  and  materials,  there  is  clearly  considerable  scope  for 
future  work  in  this  field.  In  any  case,  it  appears  that  for  less  than  10%  reinforcement 
loading,  modeling  of  nonreinforced  composites  may  be  used  to  satisfactorily  predict  lateral 
nxxle  behavior. 


6.2 


i  as  a  Function  of  Polymer  Elastic  Properties 


The  dependence  of  parasitic  mode  behavior  on  temperature  has  also  been 
investigated.  As  temperature  is  increased,  the  polymer  passes  through  its  “glass  transition” 
temperature.  At  temperatures  above  this  point  the  polymer  becomes  soft  and  highly 
attenuative.  Generally,  polymers  which  are  relatively  soft  at  room  temperature  have  lower 
glass  transition  temperature.  The  impedance  characteristics  of  the  baseline  sample  (#16822) 
was  tested  from  -SOP  to  -f50°C  at  10°  intervals.  These  impedance  spectra  are  illustrated  in 
Fig.  56.  It  is  readily  evident  that  the  glass  transition  teiiq)erature  of  the  polymer  used  in  this 
sample  is  approximately  +30°C.  This  is  relatively  low  for  an  epoxy  polymer.  Many 
applicaticms  would  require  that  the  glass  transititxt  temperature  be  well  outside  the  qjCTating 
temperature  range  of  the  practical  device.  As  the  ten^rature  increases,  the  frequencies  of 
the  resonances  reduce  as  the  polymer  becomes  more  compliant.  Careful  examination 
reveals  that  the  decrease  in  frequency  in  the  resonant  frequencies  is  greater  for  the  lateral 
modes  than  for  the  thickness  mode.  This  is  probably  due  to  the  fact  that  the  polymer’s 
shear  wave  velocity  decreases  more  rapidly  than  its  longitudinal  wave  velocity.  This 
implies  that  the  Poisson's  ratio  of  polymer  increases  as  polymer  softens.  This  is  a 
satisfying  result  since  it  is  well  known  that  compliant  polymers  usually  possess  higher 
PcHSSon’s  ratios  than  less  compliant  polymers.  The  major  feature  observed  as  the  polymer 
passes  through  the  glass  transition  temperature  is  that  the  polymer  becomes  so  attenuative 
that  it  can  no  longer  support  lateral  resonant  activity.  Although  this  may  appear  to  be  a 
desirable  characteristic,  it  is  also  probable  that  thickness  energy  is  not  being  coupled 
efficiently  between  the  active  thickness  vibrating  ceramic  pillars  and  the  neighboring 
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17493. 

A:  IZl  B:  $  o  MKR  259  750.000  Hz 

A  MAX  10.00  KQ  MAG  436.467  Q 

B  MAX  80.00  deg  PHASE  -17.0861  deg 


A  MIN  200.0  0  START  100  000.000  Hz 

B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Figure  55 

Inyedance  characteristic  of  sample  #17493  (Ceramaform  reinforced  sample) 
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8  o  mKR  307  569.250  Hz 

A  MAX  1.000  f1C  MAG  1.53219  KO 

B  MAX  60.00  deg  PHASE  23.6250  deg 


-SO  DEG  C  _ 

PM  16S22B 

A:  IZI  B:  «  o  MKR  305  069.500  Hz 


A  MIN  250.0  S  START  100.000  Hz 
B/DIV  20.00  deg  STOP  1  000  000.000  Hz 
-40  DEG  C  _ 


Figure  56 

Impedance  characteristic  of  sanq}le  #16822  (Baseline)  measured  at  various 
ten^eratures  (10°C  intervals) 
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PN  leeaaB 

a:  IZI  B:  o  MKR  300  070.000  Hz 

A  MAX  1.000  Mfl  MAG  1.64458  Kft 

B  MAX  60.00  dea  PHASE  29.8430  deg 


A  MIN  250.0 

Q 

START 

100.000  Hz 

B/OIV  20.00 

deg 

STOP 

1  000  000.000  Hz 

-30  DEG  C 

PN  168226 

A:  IZI  B:  0 

o  MKR 

292  570.750  Hz 

A  MAX  1.000 

Mfi 

MAG 

1 . 44646  Kfi 

B  MAX  60.00 

deg 

PHASE 

34.3040  deg 

A  MIN  250.0  S  START  100.000  Hz 
B/DIV  20.00  deg  STOP  1  000  000.000  Hz 
-20  DEG  C 


Figure  56 
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aB2  571.750  Hz 
1.41460  Kfi 
38.9210  d»g 


PM  1632BB' 


A:  IZJ  B:  ^  o  MKR 

A  MAX  1.000  Mfi  MAG 

B  MAX  60.00  d*g  PHASE 


-10  DEG  C  _ 


PN  1G622B 

a:  IZi  B:  B  o  (IKR  280  072.000  Hz 

A  MAX  1.000  Mfi  riAG  1.4515B  Kfi 


0  DEG  C 


Figure  56 
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PN  isaasB 

ft:  IZI  B:  e  O  MKR  272  572.750  Hz 

P  MAX  1 . 000  MQ  MAG  1 . 36966  KQ 


B  MAX  60.00  deg  PHASE 

P  *  ‘  ^  !  • 


39.7636 


I  /ji  !  =  _  *  ^ 

■■■''■Tnr“i  I"  “1  ! 


"I - p- 


I  I 


A  MIN  250.0  fi  START  100.000  Hz 
B/OIV  20.00  deg  STOP  1  000  000.000  Hz 
+10  DEG  C_ 


PN  16S22B 

A:  IZI  Bi  0  o  MKR  262  573.750  Hz 

A  MAX  1.000  M«  MAG  1.66211  Kfl 


A  MIN  250.0  12  START  100.000  Hz 
B/DIV  20.00  deg  STOP  1  000  000.000  Hz 
+20  DEG  C 


Figure  56 
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347  575.350  Hz 
i  .  37966  KSl 
34.4306  dvg 


PN  16833E 

A:  iZI  B:  o  HKR 

A  MAX  1 . 000  MS  MAG 


A  MIN  350.0  «  START  100.000  Hz 
6/DIV  30.00  d«g  STOP  1  000  000.000  Hz 
+3Q  DEG  C 


PN  168228 

a:  121  B:  ^  o  MKR  *  235  076.500  Hz 

A  MAX  1.000  MR  MAG  984.156  » 


A  MIN  250.0  fl  START  100.000  Hz 
B/OIV  20.00  deg  STOP  1  000  000.000  Hz 
+40  DEG  C  _ 


Figure  56 
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PN  1 fiRPPR 

n:  IZI  B:  a  o  MKR  247  575.250  Hz 

A  MAX  1  . 000  MiJ  MAG  1 . 40679  Kfl 


B  MAX  6Q.00  deg  PHASE  19.6273  deg 


A  MIN  250.0  START  100.000  Hz 
B/DIV  20.00  deg  STOP  1  000  000.000  Hz 
+50  DEG  C_ 


Figure  56 
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polymer.  Efficient  composite  transducer  operation  relies  upon  efficient  coupling  of 
thickness  polarized  vibration  between  the  pillars  and  the  surrounding  polymer  matrix. 

6.3  Performance  of  “Enhanced”  Comix>sites 

The  motivation  for  designing  1:3  composites  with  aligned  vertical  compressible 
inclusions  has  been  discussed  at  length  in  Section  S.  FMI  manufactured  three  enhanced 
cmnposites.  Each  composite  was  split  into  quadrants,  each  possessing  a  different  hole 
getmietry.  These  gecmietries  are  illustrated  in  Fig.  57.  Sample  #16978  (Fig.  57)  is  similar 
to  the  baseline  design  (#16822)  but  possesses  0.020"  holes.  Sample  #16979  has  the  same 
patterrL  but  with  0.040"  holes.  The  impedance  sp«;tra  of  the  quadrants  of  these  composites 
are  iUustrated  in  Figs.  58  and  59  respectively.  It  is  clear  that  the  0.020"  holes  have  had  little 
impact  on  the  resonant  characteristics  present  in  the  undrilled  sample.  However,  when  the 
results  for  #16979  are  considered,  it  is  evident  that  the  larger  holes  have  reduced  the 
strength  of  the  higher  frequency,  diagonal,  parasitic  lateral  mode.  It  was  expected  that  the 
presence  of  holes  would  severely  disrupt  the  lateral  resonant  characteristic.  The  phenomena 
which  give  rise  to  these  results  are  not  fully  understood  and  further  investigation  would  be 
valuable.  In  any  case,  samples  #16978  and  #16979  were  designed  for  use  at  low 
frequencies  where  parasitic  resonant  behavior  is  immaterial.  The  devices  were  sent  to  NRL 
USRD  for  hydrostatic  sensitivity  testing.  No  results  are  available  at  this  time. 

A  third  enhanced  composite,  with  a  ceramic  volume  fraction  of  30%  was 
manufactured  by  FMI.  This  was  tested  at  Stanford  and  the  results  obtained  presented  in 
Section  5  (Table  I).  A  modest  improvement  in  thickness  rtKxle  transducer  sensitivity  in  both 
transmission  and  reception  was  observed  when  the  quadrants  possessing  compressible 
inclusions  were  used. 

The  work  on  enhanced  composites  is  at  an  early  stage.  Initial  theoretical  and 
experimental  work  suggest  the  promise  of  improved  efficiency.  This  is  most  evident  when 
a  low  elastic  modulus,  high  Poisson's  ratio  polymer  is  used.  Additionally,  greatest  benefits 
are  expected  when  the  transducer  is  used  in  the  low  frequency  hydrophone  mode.  The 
economics  of  rruuiufacture  are  yet  to  be  addressed.  However,  there  are  some  techniques 
currently  being  developed  which  hdd  promise  for  low  cost  manufacture. 
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PZT  COMPOSITE  DESIGN 
16978 


EACH  QUADRANT  TO  H/WE  8X8  PZT  ROD  ARRAY 


Figure  57 

Configuration  of  holes  drilled  in  the  quadrants  of  sample  #16978 
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ft  MIN 
B/DIV 


508.  Q 

20.00 


START 

STOP 


000  000.000  Hz 


PN  1697S~S 

A:  IZI  B:  ^  o  MKR  257  500.000  Hz 

A  MAX  50.88  Kfi  MAG  4.00004  Kit 

B  MAX  00.00  dtfl  PHASE  38.1890  deg 


A  MIN  500.0  fi  START  100  000.000  Hz 
6/D  IV  20.00  deg  STOP  1  000  000.000  Hz 


Figure  58 

Impedance  characteristic  of  the  quadrants  of  sample  #16978  (0.020"  holes) 
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PN  15578-5 

A:  IZl  B:  o  mKR  262  900.000  Hz 

A  MAX  50.00  KJI  MAG  4.44958  KSJ 

B  MAX  80.00  deg  PHASE  41.3082  deg 


A  MIN  503.0  ri  START  100  800.380  Hz 

B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


PN  16578-4 

A;  121  B:  e  o  MKR  259  750.000  Hz 

A  MAX  50.00  Kft  MAG  3.92204  Kfl 

B  MAX  80.00  deg  PHASE  44.6725  deg 


A  MIN  500.0  2  START  100  000.000  Hz 
B/DIv  20.00  deg  STOP  1  000  000.000  Hz 


Figure  58 
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PN  16579-1 

A*.  IZl  B:  o  MKR  259  750-000  Hz 

A  MAX  50.00  Kfi  MAG  5.43732  KK 

B  MAX  80.00  deg  PHASE  36.4047  deg 


AMIN  500 . 0 
B/DiV  50.00 


S  START 
deg  STOP 


100  000.080  Hz 
1  000  000.000  Hz 


PN  16979-2 
A-.  121  B:  $ 

A  MAX  50 . 00 


KR 


o  MKR 
MAG 


253  000-000  Hz 
4.74197  KR 


Figure  59 

Irapedancc  characteristic  of  the  quiulrants  of  sample  #16979  (0.040"  holes) 
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PN  16975-5 

A:  IZI  B:  ^  c>  MKR  390  000,000  Hz 

A  MAX  50,00  KJJ  MAG  6.64733  KJJ 

B  MAX  60.00  deg  PHASE  30.7866  deg 


A  MIN  500.0  ^  START  100  000.000  Hz 

B/DIV  60.00  deg  STOP  1  000  000.000  Hr 


PN  16979-4 

A:  IZI  B:  o  mKR  268  000.000  Hz 

A  MAX  50.00  KQ  MAG  4.46121  Xft 


A  MIN  500.0  51  START  100  000.000  Hz 

B/DIV  20.00  deg  STOP  1  000  000.000  Hz 


Figure  59 


6.4  Cgnclusions 


During  the  course  of  this  collaborative  program  it  has  become  clear  that  the  concept 
of  using  distributed  period,  and  specialized  geometry,  designs  is  successful  as  a  means  of 
reducing  parasitic  resonant  behavior.  However,  it  would  appear  that  it  is  still  desirable  to 
use  the  finest  scale  structures  available  within  bounds  set  by  economic  and  manufacturing 
factors.  The  temperature  tests  have  yielded  an  improved  understanding  of  the  different 
modes  and  their  dependence  on  polymer  elastic  modulus.  In  general  it  is  desirable  to 
employ  a  high  modulus  polymer  with  a  high  glass  transition  temperature.  A  limited  amount 
of  work  with  fiber  reinforcement  has  proven  inconclusive  at  this  stage.  It  is  foreseeable  that 
the  rewards  of  improved  mechanical  stability,  enhanced  sensitivity  and  reduced  resonant 
activity  suggest  plentiful  scope  for  fruitful  woric  in  future.  Finally,  the  work  on  “enhanced” 
composites  has  proven  very  interesting  and  produced  encouraging  initial  results.  There  is 
an  obvious  motivation  for  extending  the  currently  limited  understanding  of  the  actual 
performance  enhancement  that  can  be  achieved  at  an  economical  cost.  There  is  also  a 
requirement  for  identifying  new  applications  for  the  large  scale  composites  that  have  been 
studied.  These  include  the  non  destructive  examination  of  lossy  media  (concrete,  wood, 
GRP  etc.),  high  frequency  sonar  (civil  and  military  applications)  and  as  actuators  in 
responsive  structures,  for  vibration  control. 

7.  MODELING  OF  TWO-DIMENSIONAL  DISTRIBUTED  PERIOD 
COMPOSITES 

In  a  standard  circular  disk  thickness  transducer,  consisting  of  an  electroded 
piezoceramic  disk,  the  main  resonance  is  a  A/2  thickness  mode  excited  directly  by  the  z- 
directed  E  field  due  to  the  vobage  applied  to  the  electrodes.  At  the  same  time,  lateral 
(radial,  in  this  case)  modes  are  excited  through  the  thickness  vibrational  motion,  by  the 
Poisson  effect  at  the  free  edge  of  the  disk.  These  lateral  modes  are  of  extensional  type, 
driven  by  the  radial  bulging  in-and-out  of  the  outer  circumference  of  the  disk.  Lateral 
modes  are  not  efficiently  excited  by  the  radially  uniform  z-directed  E  tield  driven  by  the 
electrodes,  because  of  spatial  averaging  over  the  radial  standing  waves  that  constitute  the 
lateral  modes.  In  this  type  of  transducer  the  lateral  modes  may  be  called  dimensional 
modes,  because  they  have  frequencies  determined  by  the  radial  dimension  of  the 
transducer. 

In  a  1-3  composite  transducer  the  piezoelectric  is  distributed  in  piezoceramic 
elements  urans^  in  either  a  strictly  periodic  or  a  distributed  period  lattice  array  that  is 
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embedded  in  a  polymer.  The  piezoelements  may  be  arranged  periodically  in  a  simple 
square  array,  or  in  rectangular,  triangular,  etc.,  periodic  arrays.  Other  possible 
configurations  have  additional  piezoelements  in  each  unit  cell  of  the  lattice. 

The  composite  structures  described  above  improve  the  characteristics  of  the 
thickness  mode  resonance.  But  they  also  alter  the  lateral  mode  spectrum.  First,  these 
structures  present  stopbands  to  the  laterally  propagating  waves,  preventing  the  formation  of 
standing  wave  resonances  at  frequencies  within  a  stopband.  Second,  at  frequencies  outside 
the  stt^bands  the  standing  lateral  waves  are  coupled  to  the  electrodes  only  at  discrete  spatial 
positions.  This  means  that  the  spatial  cancellation  of  lateral  mode  coupling  that  occurs  in  a 
homogeneous  ceramic  transducer  is  lifted  in  a  strictly  periodic  composite  transducer.  In 
particular,  at  the  stopband  edges  the  ceramic  elements  all  vibrate  in  phase,  or  alternatively, 
in  antiphase.  The  in-phase  bandedge  standing  waves  are  strongly  coupled  to  the 
electrodes,  while  the  antiphase  bandedge  standing  waves  are  totally  decoupled. 

To  summarize,  then,  in  a  1-3  periodic  composite  some  of  the  lateral  dimensional 
resonances  are  suppressed  by  the  stopband  structure  of  the  lateral  wave  dispersion  curves. 
On  the  other  hand,  certain  lateral  dimensional  resonances  at  the  stopband  edges  (the  in- 
phase  bandedge  standing  waves)  experience  greatly  enhanced  coupling  to  the  electrodes. 
These  strongly  coupled  lateral  resonances  degrade  the  performance  of  the  transducer  if  they 
are  located  too  close  to  the  thickness  resonance  frequency. 

The  conventional  approach  to  the  problem  of  bandedge  resonances  is  to  reduce  the 
period  of  the  periodic  array,  shifting  the  resonances  to  frequencies  much  higher  than  the 
thickness  resonances.  Another  solution  is  to  introduce  a  quasirandom  distribution  in  the 
spatial  periodicity  of  the  ceramic  cylinder  elements.  A  simple  one-dimensional  analytic 
OKKlel  (based  on  transmission  line  theory)  has  shown  that  a  large  one-dimensional  array 
having  a  distribution  of  periods  exhibits  a  large  reflective  attenuation  (Section  2.3), 
suggesting  that  such  an  array  would  effectively  suppress  all  lateral  resonances. 
Experiments  on  two-dimensional  arrays  with  a  Gaussian  distribution  of  the  unit  cell 
dimensions  have  demonstrated  that  the  bandedge  resonances  can  be  effectively  controlled 
using  this  method — which  tends  to  suppress  the  lateral  resonances  by  virtue  of  the 
increased  lateral  wave  attenuation.  It  also  tends  to  reduce  piezoelectric  excitation  of  these 
modes  by  randomizing  the  phases  of  the  coupling  to  the  electrodes  through  the  individual 
ceramic  cylinders. 
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Up  to  now  the  design  of  two-dimensional  distributed  period  1-3  composite 
transducers  has  been  almost  entirely  empirical.  Floquet  Theory  is  applicable  only  to 
periodic  arrays;  and  application  of  the  finite  element  method  to  a  large  two-dimensional 
array  is  impractical,  as  it  would  require  computation,  unless  the  domain  of  computation 
reduced  to  a  single  periodically  repeated  unit  cell.  For  improved  modeling  of  distributed 
periodic  arrays  there  is  need  of  an  alternative  approach.  There  appears  to  be  no  analytic 
method  analogous  to  transmission  line  theory  for  one-dimcnsional  distributed  period 
composites.  But  another  approach  having  advantages  and  limitations  similar  to  those  of 
Floquet  Theory  is  based  on  a  scattering  analysis  formulation  using  the  Reciprocity  Relation 
for  elastic  fields. 

The  basic  idea  is  to  use  the  elastodynamic  reciprocity  relation  to  calculate  the 
reflection  and  transmission  through  a  slab  of  composite  material.  In  particular,  this  permits 
calculation  of  the  transmission  loss  due  to  scattering  by  the  ceramic  elements  of  the 
composite.  The  transmission  line  model  referred  to  above  can  be  used  to  calculate 
transmission  loss  for  one-dimensional  composites,  but  not  for  the  two-dimensional  case. 
Using  the  reciprocity  approach,  both  cases  may  be  treated. 

Scattering  theory  based  on  the  reciprocity  relation  calculates  the  difference  in 
scattering  for  two  different  geometries.  In  the  problem  at  hand  the  two  geometries  are; 

(i)  a  uniform  polymer  medium  (unpiimed  quantities) 

(ii)  a  composite  slab  embedded  in  polymer  (primed  quantities). 

The  change  in  scattering  can  then  be  expressed  as 

A5  =  -^  /  {(p'  -  P)v'  •  V  -  T':  (s'  -  s):  T}  dV 

^  */ ceramic  '•  ^ 


where  the  primes  on  the  density  p  and  the  compliance  s  refers  to  the  ceramic  properties. 
To  calculate  reflection,  the  primed  velocity  and  stress  fields  are  for  wave  incidence  on  the 
same  side  of  the  slab.  For  a  transmission  calculation  the  unprimed  fields  are  for  wave 
incidence  on  the  left  side  of  the  slab,  while  the  primed  fields  correspond  to  incidence  on  the 
right.  Since  geometry  (i)  has  zero  reflection  and  complete  transmission,  the  reflection  and 
transmission  for  geometry  (ii)  can  both  be  evaluated  from  AS. 
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To  illustrate,  the  reflection  of  shear  waves  from  a  slab  of  one-dimensional  (2-2) 
con^)osite  can  be  shown,  for  the  low  volume  fraction  case,  to  be 


R  = 


44  L 


P'  -P  + 


pC. 


44 


^4  -  C. 


44 


,  A’-l 

E 

■'  n=0 


UK 


(29) 


for  a  slab  consisting  of  N  ceramic  elements  at  positions  Zn,  with  thicknesses  Wn. 
Printed  material  constants  are  for  the  ceramic,  and  unprimed  for  the  polymer. 

As  stated,  the  formula  applies  to  a  general  distributed  period  structure  with 
nonuniform  ceramic  elements.  For  the  strictly  periodic  case 

Wn  =  to 


In  the  second  case,  scattering  from  the  individual  ceramic  element,  detailed  by  the 
summation,  adds  coherently  to  produce  stopbands  and  bandedge  resonances.  For  a 
distributed  period  structure  this  coherence  is  broken,  eliminating  the  stopbands.  At  the 
same  time,  a  broadband  transmission  loss  results  from  the  incoherent  scattering, 
ccKiesponding  to  the  transmission  coefficient 

ir|  =  y/\  -  jiip 


The  advantage  of  the  reciprocity  approach  is  that  it  can  be  easily  extended  to  two- 
dimensional  (1-3)  composites.  In  this  case  the  calculation  is  simplified  by  taking  the  two 
getxnetries  in  the  reciprocity  relation  to  be: 

(i)  a  strictly  periodic  1-3  composite  slab 

(ii)  a  distributed  period  1-3  con^site  slab. 

Transmissitm  loss  for  case  (i)  can  be  calculated  by  using  Floquet  Theory;  and  from  this 
result  the  transmission  loss  for  case  (ii)  can  be  calculated  using  AS  obtained  from 
reciprocity. 
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A  brief  introduction  to  this  work  was  given  at  the  ONR  Review  of  piezoelectric  and 
electrostrictive  materials  to  transducer  materials.  State  College,  PA,  April  1, 1992.  It  will 
be  elaborated  in  a  projected  review  paper  on  1-3  composites,  theory  and  practice. 

8.  CONCLUSIONS 

Theoretical  work,  supported  by  experimental  results,  has  proven  that  parasitic 
lateral  mode  behavior  in  1:3  composite  transducers  may  be  suppressed  by  using  a  statistical 
distribution  of  inter-pillar  periods.  Alternatively,  a  similar  effect  can  be  obtained  by  using 
double  or  triple  p>eriods.  Factors  crucial  to  the  successful  manufacture  of  0:3  transducers 
have  been  identified.  These  include  the  use  of  clean,  rounded  particles,  possessing  a 
narrow  size  distribution  compressed  under  high  pressure.  A  simple  model,  taking  account 
of  all  the  relevant  factors,  has  been  developed.  This  represents  a  signiOcant  improvement 
over  previously  used  models.  Nevertheless,  the  mechanisms  of  ferroelectricity  near  the 
surfaces  of  piezoelectric  ceramic  particles  present  complex  issues  requiring  further  study. 

The  concept  of  using  multiple,  independently  excited,  layer  transducers  has  been  proposed 
and  tested.  Simple,  yet  comprehensive,  modeling  techniques  have  been  developed.  It  has 
been  shown  theoretically  and  experimentally  that  these  transducers  are  capable  of  hitherto 
unheard  bandwidth  performance.  Recently,  1:3  transducers  with  compressible  gas  filled 
inclusions  have  been  investigated.  The  inclusions  absorb  lateral  strains  which  would 
otherwise  couple  in  an  undesirable  manner  to  thickness  polarized  charge  in  the  ceramic 
pillars.  There  is  reason  to  believe  that  a  significant  improvement  in  hydrophone 
performance  may  be  achieved  using  this  technique. 
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IS.  APPENDICES 


Appendix  I;  Propgnics  of  PZT  5A.  Epoxy  and  Polyurethane 


PZT  5  A 

Epoxy 

Polyurethane 

ciiE(109Nm-2) 

121 

9.63 

9.01 

ci2^  (109  Nm-2) 

75.4 

5.2 

7.99 

C33E  (109  Nm-2) 

111 

9.63 

9.01 

c44E(io9Nm-2) 

21.1 

2.22 

0.51 

p  (kgm-3) 

7750 

1140 

1300 

e3i  (Cm-2) 

-5.4 

0 

0 

e33  (Cm-2) 

12.3 

0 

0 

ei5  (Cm-2) 

12.3 

0 

0 

eil'^/60 

1730 

4 

4 

e33^/eo 

1700 

4 

4 

Appendix  U:  Design  parameters  of  the  double  layer  transducers 


TRANSDUCER 

A 

B 

C 

D 

Diameter  (mm) 

12 

12 

18 

18 

Layer  thickness  (mm) 

1.0 

1.0 

0.5 

0.5 

Active  imterial 

Lead  Titanate 

Lead  Titanate 

Composite 

Composite 

Backing  (MRayl) 

10.0 

20.0 

10.0 

16.4 

NOTE:  The  composite  material  was  50%  PZT-5H  type  ceramic.  Transducers  A  and  C  were 
fabricated.  The  backings  were  composed  of  tungsten  loaded  epoxy.  The  layers  were 
bonded  via  a  fine  epoxy  layer.  No  matching  layer,  other  than  a  thin,  protective  coat  of 
paint,  was  used. 
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Parameter 

Transducer  I 

Transducer  II 

Ceramic  volume  firactitxi  (%) 

20 

30 

Thickness  (mm) 

2.5 

3.8 

Pillar  width  (mm) 

0.75 

1.00 

Hole  diameter  (mm) 

0.5 

0.5 

Protective  layer  thickness  (mm) 

0.3 

1.0 

Center  frequency  (MHz) 

0.5 

0.35 

Notes:  The  ceramic  is  PZT-5H.  Transducer  A  employed  Hysol  RE2039-HD3561  epoxy. 
Transducer  B  was  manufactured  by  FMI,  Maine,  and  used  a  proprietary  epoxy. 
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